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ABSTRACT 

We measure the projected 2-point correlation function of galaxies in the 180 deg^ equato¬ 
rial regions of the GAMA II survey, for four d i fferen t redshift slices between z = 0.0 and 
z = 0.5. To do this we further develop the lColel (1201 ll) method of producing suitable random 
catalogues for the calculation of correlation functions. We find that more r-band luminous, 
more massive and redder galaxies are more clustered. We also find that red galaxies have 
stronger clustering on scales less than ~ 3/i“^Mpc. We compare to two different versions 
of the GALFORM galaxy formation model, Lacey et al (in prep.) and iGonzalez-Perez et alJ 
(|2014|) . and find that the models reproduce the trend of stronger clustering for more massive 
galaxies. However, the models under predict the clustering of blue galaxies, can incorrectly 
predict the correlation function on small scales and under predict the clustering in our sample 
of galaxies with ~ 3L*. We suggest possible avenues to explore to improve these cluster¬ 
ing predictions. The measurements presented in this paper can be used to test other galaxy 
formation models, and we make the measurements available online to facilitate this. 

Key words: large-scale structure of Universe, galaxies: evolution, galaxies: formation 


1 INTRODUCTION 


The two-point auto-correlation function is a widely used statistical 
description of the spatial distribution of galaxies. In the standard 
model of cosmology, the shape of this function is set both by the 
mass of dark matter haloes in which a particular galaxy sample re¬ 
sides, and baryonic processes which can change the spatial distribu¬ 
tion of galaxies on smaller scales. Its dependence on galaxy proper¬ 
ties is well established, at low redshifts, by large area spectrosc opic 
surveys such as SDSS jYork et aD l200(T Istrauss et ai.ll2002l) and 
2 Degree Field Galaxy Redshift Survey (2dFGRS; Colless et~d] 
I2OOII1 . The amplitude of the auto-correlation function of galax¬ 
ies is seen to be strongly dependent on luminosity, stellar mass 
and colour. At low redshifts, brighter, redder and more massive 
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galaxies have been observed to be more strongly clustered (e.g . 


NorbergetalJl2001.2002l: 

Li et alJl2006l:IZehavi et al.ll200ll200l 

201lLIChristodoulou et all 

2012h. 


In the higher redshift Universe, small (< 1 deg^) but deep 
spectroscopic surveys have also measured the clustering of galax¬ 
ies. The DEEP2 survey has been used to demonstrate that the colour 
dependence of galaxy clustering is already in place at z ~ 1, 
whilst within a red or blue sample of galaxies, clusterin g is insen¬ 
sitive to luminosity over the range 20.2 < Mb < 21.8 l lCoil et alj 
l2008h . Compared with lower redshift SDSS data, the DEEP2 mea¬ 
surements of the clustering of brighter and more massive galax¬ 
ies has a larger amplitude than expected from scaling the low red¬ 
shift measurements using linear theory dCoil et al.ll2008l : [ui et alj 
I2OI2I) . This can be interp reted as evidence of significant bias evo¬ 
lution for these galaxies dCoil et alJl2Q08l : lLi et alj|2012ll . Another 
example of a small area, deep spectroscopic survey is the VIMOS- 
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VLT Deep Survey (VVDS), which found a sharp increase in the 
amplitude of galaxy clustering around the characteristic magni¬ 
tude of the sample’s luminosity function dPollo et al.l2006h . VVDS 
also found, in agreement with DEEP2, that the clustering of mas¬ 
sive galaxies, with stellar mass > at 2 ; ~ 1 can 

only he reconciled wi th lower redshift me asurements if their bias 
evolves significantly dMeneux et alJ[20Q3) . Results from another 
small area, deep survey, zCOSMOS, sho w that this is true for g alax- 
ies more massive than > dMeneux et alJIlOO^ . The 

VIPERS survey has also found that more luminous and more mas¬ 
sive galaxies are more clustered at redshi fts of 0.5 < 2 < 1.1 
than their fainter, le ss massive counterparts dde la Torre et alj2013l: 
iMarulli et alj 2013h. More recently, the 9 deg^ PRIMUS survey 

1 Coil et alj|201 ll : Cool et a02013h has found red galaxies are more 
clustered than blue galaxies, and the relationship between cluster¬ 
ing amplitude and sample luminosity is in place in the re dshift 
slices; 0.2 < 2 < 0.5 and 0.5 < 2 < 1.0 dSkibba et alj2014h . 

In addition to spectroscopic surveys, surveys relying on pho¬ 
tometric redshifts have also explored the clustering of galaxies. 
Galaxies with red colours have been shown to have steeper cor¬ 
relation functions th an bluer galaxies at 2 ~ 0.5 in the 0.78 deg^ 
COMBO-17 survey dPhleps et~^l2006h . and the 1.5 deg^ Ultra- 
VISTA survey demonstrated that m ore massive galaxies are more 
clustered between 0.5 < 2 < 2.5 dMcCracken et al.ll2014h . Pho¬ 
tometric surveys have also studied the evolution of ga laxy clus¬ 
tering. In the 6.96 deg^ Bootes field iBrown et alj l l2008h found, at 
0.2 < 2 < 1.0, that red galaxies fainter than Mb =20 — 5 logj^Q h 
showed no luminosity dependence of the amplitude of their cluster¬ 
ing, whilst the converse is true for red galaxies brighter than this. 
Studies in the Bootes field also demonstrated luminous red galaxies 
display little evolution in the amplitude of their clustering be tween 

2 = 0.5 and 2 = 0.9 dWhite et al.ll2007l : [Brown et alJl2008h . This 
observation, which suggests the clustering of these galaxies evolves 
slower than the underlying dark matter distribution, can be ex¬ 
plained by the rem oval of highly-bias e d satellite galaxies by merg¬ 
ing or disruption dWhite et alj l2007l : iBrown et'Sl l2008h . Whilst 
bright, red galaxies display little evolution in clustering strength 
with redshift, fainter galaxies (—22 < Mb — hlog^^ < —19) 
have been shown to have a decreasing clustering amplitude between 
2 = 0.4 and 2 = 1.2 in the CFHTLS photometric redshift survey 
dMcCracken et iDIIooi). 

The photometric surveys at high redshift are still fairly small 
(< 10 deg.^). A possible exception to the trade-off betw een large 
area and high redshift is the work of IGuo et alj d2014l) ; here the 
authors used the CMASS sample of SDSS to show that clustering 
was stronger for brighter and redder galaxies at 2 ~ 0.5. However, 
the complicated selection of the CMASS sample limits their work 
to a very narrow luminosity and redshift range. 

The GAMA spectroscopic survey offers a new window onto 
the clustering of galaxies and its evolution with redshift. It has 
a larger area (the equatorial regions we use total 180 deg^) than 
the deep surveys but has spectra of much fainter galaxies (2 mag¬ 
nitudes) than large-area surveys like SDSS. As such it comple¬ 
ments both types of survey, by enabling detailed clustering mea¬ 
surements at an intermediate epoch. Whilst our redshift range over¬ 
laps with the PRIMUS low-redshift sample, our larger area allows 
us to split the data into finer redshift bins and consider brighter, 
rarer objects. In this paper we study the projected two-point corre¬ 
lation function (2PCF) of galaxies as a function of their luminosity, 
stellar mass and colour, in four different redshift bins, in order to 
fully appreciate any redshift evolution. To aid in the physical inter¬ 
pretation of our results we compare them to different versions of 


the semi-analytic galaxy for mation model GALFORM, specifically 
iGonzalez-Perez et al.l l l2014h (hereafter G14 model) and Lacey et al 
(in prep.) (hereafter L14 model). We quantitatively compare the rel¬ 
ative merits of the different versions of the GALFORM model and 
discuss the impact of this on our understanding of galaxy forma¬ 
tion. To see a complementary approach, we refer the reader to Pala- 
mara et al. (in prep.), who fit halo occup ation distribution models 
(HOD models; see e.g. lZheng et alJlloO^ to GAMA. 

This paper is organised as follows. Section introduces the 
GAMA data and the galaxy formation model we use, along with 
details of how we calculate luminosity, mass and rest-frame colour. 
Section presents our method of generating a random catalogue 
and measuring clustering. In Section |4] we present our results, be¬ 
fore discussing them and concluding in Section |5] We present our 
results using units of a fiducial ACDM cosmology with flm = 
0.25, Da = 0.75 and H = lOOh km s'^ Mpc"^ 


2 DATA AND MODELS 

2.1 The Galaxy and Mass Assembly survey 

The GAMA survey is a spectroscopic and multi-wavelength sur- 
vey of galaxies carried out on the Anglo-Australian telescope 
dPriver et al.ll201 ll : iLiske et alJl^lSll . In this work we utilise the 
main r-band limited data from the GAMA II equatorial regions, 
which consists of a highly complete (> 98%) spectral catalog ue 
of galaxies selected from the SDSS DR7 dAbazaiian et alj|2009f l to 
have Fpetro < 19.8. The older GAMA I survey had a shallower 
limit of r < 19.4 in two of the regions. Target fields are repeat- 
edly observed in a wa y that removes biases against close pairs 
dRobotham et alj|2010h . avoiding such biases is ideal for cluster¬ 
ing measurements. Star/galaxy separation is based on the difference 
between an object’s model and PSF magnitudes in SDSS DR7 data 
and, where UKIDSS photometry is available, the object’s optical 
a nd infrared co l ours. Further details of the GAMA survey are given 
in fBaldrv et aLl dTOf^ . lRobotham et al] d2010h , [Driver et al] d201lh 
and lLiske et al.l d2015h . The data are split over three 12 x 5 deg^ 
fields centered at 9'* (G09), 12'‘ (G12) and 14.5'* (G15) R.A. and 
approximately 5 = 0 degrees declination. 


2.1.1 Redshifts 

Redshifts for GAMA objects were measur ed automatic a lly, us ¬ 
ing the so f tware AUTOZ as described in iBaldrv et"^ d2014b . 
iLiske et H] d2015h find 0 .2% of the sample are expected to have 
an incorrect redshift. The me dian velocity unce rtainty of the mea¬ 
sured redshifts is 27 kms“'^ dLiske et al.ll2015h . The redshifts are 
taken from a table called DistancesFramesv12 in a GAMA 
Data Management Unit (DMU). These reds hifts have been cor - 
rected for the local flow using the model of iTonrv et al.l ( l2000ll . 
smoothly tapered to the CMB restframe for 2 ^ 0.03. 


2.1.2 Quality cuts 

In addition to the redshift quality cut, nQ ^ 3, we also 

only consider galaxies in regions with completeness greater than 
80% using the GAMA angular completeness mask dPriver et alj 
I 2 OIIL ILiske et ^l2015h . We additionally only select objects with 
VIS_CLASS= 0, VIS_CLASS= 1 or VIS_CLASS= 255, which 
removes objects which upon visual inspection do not show any 













































































Galaxy and Mass Assembly (GAMA): Projected Galaxy Clustering 3 



light or appear to be part of another galaxy 


2.1.3 Magnitudes 

GAMA combines data from a large number of ground and space 
based telescopes and so has a very wide wavelength range, from 
X-ray to radio. In this work we use optical photometry from SDSS 
DR7 imaging data. To define lumino sity samples we use SDSS Pet¬ 
rosian magnitudes ( |Petrosianlll976h . as the GAMA selection used 
Petrosian magnitudes. To define colours, and when estimating stel¬ 
lar mass, we use SDSS model magnitudes as these are often more 
suitable for colour terms (see the SDSS DR7 photometry web- 
pag^. Following SDSS conventions, we will label model mag¬ 
nitudes using the letter associated with the bandpass in which the 
magnitude was measured. We took these magnitudes from the table 
called TilingCatv42 in a GAMA DMU, where they are repli¬ 
cated. 

GAMA also has a set of magnitudes meas ured in aperture s 
with a matched size across the different bands ( iHill et alJl20iTh . 
these magnitudes have been shown to be superior to model mag¬ 
nitudes when calculating a gal axy’s stellar mass using spectral en¬ 
ergy distribution (SED) fitting jTavlor et al.l201 ih . Our adoption of 
model magnitudes was in order to avoid the small num ber of fail¬ 
ures w hose spatial positions have not been mapped in the lHill et alj 
( I2OIII1 catalogue. Model magnitudes have been widely used to de¬ 
fine colours in SDSS publications, and should be sufficient to sep¬ 
arate red and blue galaxies in this work. 


combined with our method of generating a random catalogue, re¬ 
moves the need to correct for den sity evolution (see Section [3TT >. 
Note that in lLovedav et al.l i2015t) a different set of parameters are 
favoured. However in Appendix IaI we show that adopting the pa- 
rameter s P — 1.45 and Q = 0.81, which we took from an earlier 
draft of iLovedav et alj j2015ll . affects the correlation function in a 
way far smaller than the errors. 

The k-corrections are derived from the GAMA DMU 
K CORR_z00v04 , which was produced usin g the m ethod set out 
in iLovedav et alJ ilOl'j) . The iLovedav et alj (1201 2h k-corrections 
are fo und by using the code KCORRECT v4 2 telanton & RoweisI 
l2007h to fit each galaxy’s u, g, r, i, and 2 -band SDSS model 
magnitudes with SED templates. For many applications the max¬ 
imum redshift at which an object fulfils the selection criteria of 
the survey, 2max, is needed. For this the k-correction as a func¬ 
tion of r edshift is requi r ed for each galaxy. To enable fast com¬ 
putation, ILovedav et alj (1201 2h fit a 4th order polynomial to the 
k-correction of each galaxy as a function of redshift. The rms dif¬ 
ference between the KCORRECT estimates of k-correction and the 


S iolynomial fits to the m is less than 0.01 magnitudes for all bands 
Lovedav et iDEoH. We further speed up the k-correction pro¬ 
cess b y using average polynomials from lMcN aught-Roberts et"^ 
(l2014h . who computed the median of the ILovedav et alj ~ 2012h k- 


correction polynomials for galaxies in seven, (g — r) rest frame 
colour bins, hereafter labelled [g — r)o. The benefit of this is a re¬ 
duction in noise introduced by fitting each galaxy with an individ¬ 
ual model dMcNaught-Roberts et al.ll2014h . Rest-fra me colours are 
still es timated from the individual k-corrections of iLovedav et alj 

( I2OI2I1 . 


2.1.4 k-corrections and evolution corrections 

To compute Petrosian absolute r-band magnitudes, Mr — 5 logj^g h 
(hereafter Mr.h), from the observed apparent Petrosian r-band 
magnitudes, m^, we apply both a k-correction, kr, and a luminos¬ 
ity evolution correction. The latter correction is necessary in order 
to compare a similar population of galaxies across time, as the lu¬ 
minosity of galaxies evolves. The relation we use is 

Mr.h = rUr- kr{z) -|- Q{z - Z^et) “ 5 logiA^Liz)) - 25, (1) 

where Q is the luminosity evolution parameter, Dl the luminosity 
distance in h~^ Mpc and 2ref is a reference redshift, for which we 
adopt 2ref = 0.0. The parameterisation of lumi nosity evolution 
we use is commonly adopted in the literature (e.g. lLin 'Zll 19991: 
ILovedav et al.ll20l3l . 

At this stage we will also introduce a parameterisation of 
density evolutio n common to luminosity function stu dies. The P 
parameter (e.g. iLin et alJ[l9^ : iLoveday et alJl2012L and refer¬ 
ences therein) parameterises the density evolution of a population 
of galaxies via 

0 *( 2 ) = ((>*(2 = 0 ) 10 '’-'^^" ( 2 ) 

where 4 >*{z) is the characteris tic number density of galaxies at red¬ 
shift 2 . ILovedav et alj j2012l) have fit evolving luminosity func¬ 
tions to a shallower version of the GAMA data, with a limit of 
Tpetro < 19.4 in G09 and G15. The y find that Q an d P are very de¬ 
generate (see also lFarrowl2013ll . In lLovedav et al.l ( l2015tl adopting 
a value of Q — 1.45 result ed in very little need for density evolu¬ 
tion (see also lFarr^l20 1 3h . We therefore adopt Q = 1.45, which 

^ http://classic.sdss.org/dr7/algorithms/photometry.html, accessed 19/8/14 


2.1.5 Stellar mass 

As well as luminosity and colour, we also want to use GAMA II 
to measure clustering as a function of stellar mass. We use the re- 
lation between col our and stellar mass found for GAMA I data by 
iTavlor et^ ( 1201 ih . namely 

logigMt(/i“^MQ) = 1 . 15 + 0 . 7 ( 5 —i)o—0.4(Mi—51ogioh) (3) 

where Mi is the rest frame i-band absolute model magnitude and 
(.<7 — ()o is the re st frame colour. This relation was found by 
iTavlor et al.l ( 1201 ih from individual estimates of the stellar mass 
of GAMA I galaxies, produced by fitting stellar po pulation syn¬ 
thesis models to the optical GAMA data (details in ITavlor et alj 
l201lh . These mass estimates should h ave a statistical Icr accuracy 
of around 0.1 dex iTavlor et alj|201 lb . The original relation used 
the GAMA matched-aperture magnitudes for the colours, whilst we 
use SDSS model magnitudes. We have tested that this does not add 
significant biases or scatter to our computed stellar masses. The re¬ 
sultant offset between our mass-to-light ratios and the ones utilizing 
the GAMA matched aperture colours with Eq. [3]is only 0.01 dex 
with an interquartile range of 0.1 dex and 10% of outliers with a dif- 
ference of > 0.14 dex. Note that the problem pTavlor et ^ i201lh 
discovered with masses using model magnitudes arose during full 
SED fitting, not using this colour relation, and were reported to be 
driven by the u-band, which has no influence on Eq.[^ 

2.2 Galaxy Formation Models & Lightcones 

2.2.1 The GALEORM model 

We compare our observations to a semi-analytic galaxy forma¬ 
tion model called GALEORM. The GALEORM model was first pre- 
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sented in ICole et al] j2000l) . but it grew from a number of pre¬ 
vious attempts to model galaxy formation. These models all as¬ 
sume galaxies form in dark matter haloes, and then use analytical 
prescriptions to app r oximate key galaxy formation processes (e.g . 

White & Ree^ll97k I White & FreSll99ll: iKauffmann etai]|l99i 

Cole et al.lll994) . The first GALFORM model has prescriptions for 
the cooling rate of gas in dark matter halos, a star formation rate 
based on available cold gas content, the attenuation of star forma¬ 
tion by supemovae feedback and the merging of satellite galaxies 
to the main galaxy l lCole et al.ll2000h . Many of these prescriptions 
have free parameters, controlling factors like the strength of super¬ 
novae feedback, that account for complicated physics not fully un¬ 
derstood or analytically described. These parameters are tuned such 
that the model fits a set of observational constraints at low redshift. 
None of the models used here were tuned to reproduce clustering 
measurements. 

The models require the history of when a dark matter halo 
formed and merged , sometimes called a halo’s ‘merger tree’ 
jLacev & Colel[l993) . For this purpose the models we present here 
use an N-body simu lation based on the Millennium Simulation 
dSpringel et al.l20 05h. but run in the more up-to-date WMAP7 cos- 
mology jKomatsu et al.ll201lT) . 


and ex tension of the iBower et al.l ( l2006h model. The Bower et"^ 
l l2006l) model added AGN feedback to the lCole et al.l ( I 2 OOOI) model 
in order to better reproduce galaxy colours, and to decrease the 
number of galaxies predicted to l ie in t he bright end of the lumi¬ 
nosity function. The lLagos et al.l ( 12012h model made the star for¬ 
mation rate in the disk of the galaxy proportional to the fraction 
of the gas that was molecular, by u sing the empirical star forma tion 
law of Blitz & RosolowsO ( I 2 OO 6 I) (see lLagos et alj201 ll . l2012l . for 
details). _ 

The update of the G14 model over iLagos et al.l ( 12012h is that it 
uses merger trees from WMAP7 cosmology. WMAP7 cosmology 
has (78 = 0 . 81, th is is smaller than in the WMAPl cosmology of 
iLagos et"^ (12012h {as = 0.9), so free parameters in the ph ysical 
prescriptions had to be retuned dGonzalez-Perez et al.ll20l4) . Tun¬ 
able parameters were adjusted in order for the predictions to still 
match the rest-frame bj and K-band luminosity functions at z = 0, 
and give reasonable evolution of the UV and K-band luminosity 
functions dGonzalez-Perez et alj|2014l) . 

The G14 model converts from star formation history and 
the IMF to a spectral energy distribution (SED) using an up- 
dated version of the stellar p opulation synthesis (SPS) model from 
[Bruzual A. & Charlotl dl99^ . 


2.2.2 Model Lightcones 

The positions and velocities of particles in the simulation are output 
at epochs spaced in expansion factor, hereafter “snapshots”. Mock 
catalogues of galaxies from GALFOR M runs on the N-body sim¬ 
ulation are produced using the code of iMerson et al.l d2013h . This 
works by interpolating between the snapshots to find the positions 
of galaxies when they enter the simulated observer’s lightcone, i.e. 
when the light from the galaxy reaches the observer. 

From the simulations 26 realisations of each model were cre¬ 
ated Iw changing the position and orientation of the virtual ob- 
serveiQ These multiple realisations allow us to improve our under¬ 
standing of sample variance (Section [3.3t . and minimise its effect 
on predictions by plotting mean measurements from the mock real¬ 
isations (Section|4). One downside is the simulation used has a lim¬ 
ited volume, 12.5 x 10^ (Mpc//i)^, as such the upper-intermediate 
and highest redshift slices will over-sample, over 26 realisations, 
the simulation by at least a factor of 3 and 8. Note this is the total 
oversampling when you consider a sum over all the realisations, an 
individual realisation is much smaller than the simulation box. Also 
mitigating this effect is the fact that the location and orientation of 
the observer is randomly assigned, in this way repeated structures 
will be observed at different redshifts and orientations, and be sam¬ 
pled by different galaxies. In the lower-intermediate redshift slice 
the volume of the 26 realisations is only 4% larger than the sim¬ 
ulation, which should result in very little oversampling. How the 
oversampling affects results will be discussed in later sections. 

The next sections describe the specifics of the two versions 
of GALFORM we consider. A summary of the relevant GALFORM 
models can be seen in Table [T] 


2.2.3 The G14 model 


The lGonzalez-Perez et alj ( 20141) (G1 4) mo del has the same phys¬ 
ical prescriptions as the Lagos et ^ ( 1201 2h model, which is itself 


^ These mock catalogues will be made available from the VirgoDB: 
http ://icc. dur. ac. uk/dat a/ 


2.2.4 The LI4 model 


The Lacey et al (in prep.) (L14) model also combines t he improve¬ 
ments to star formation rate calculations ma de in the iLagos et alj 
( l2012h model and the AGN feedback model of iBower et al J ( 2006h . 


It was also run on the new WMAP7 N-body simulation. Another 
difference of the L14 model i s its in itial mass function (IMF). The 
GI4 model uses a lKennicul^ ( Il983l) IMF for all stars formed. The 
L14 model instead adopts a top-heavy IMF i n starbursting galax- 
ies, as developed in, but less extreme than, the lBaugh et al.l ( 120050 
model. This change was motivated to bring the predictions of sub- 
mm galaxy number counts into closer agreement with observations 
dBaugh et alj2003l . 

A further difference in the L14 model which could have an im¬ 
pact on clustering predictions is the treatment of merging satellite 
galaxies. Th e models we pres ent here all use an analytic approxi¬ 
mation, from ICole et al.l ( l2000l) . to calculate how long it takes for an 
accreted satellite to merge with the central galaxy. This approxima¬ 
tion assumes satellites enter the halo on orbits ra ndomly s electe d 
from the distribution of satellite orbits given in iTormenI (Il997h. 


before the orbit decays due to dynamical friction. The Cole et al] 
( I 2 OOOI) method calculates the time for the orbit to decay using 
the formula for Chandrasekhar dyna mical friction in an isothermal 
sphere given in Lace^_&Cole (Il993h. T heL14 model, on the other 
hand, uses the I Jiang et alj ( 20oS ^ 2014h formula for the timescale. 
The formula has been empirically modified using N-body simu¬ 
lations in order t o account for the tida l stripping of the accreting 
haloes. Note tha t Campbell et alj (l2015h . motivated by the work of 
IContreras et alj ( l2013h . additionally develops the model to track the 
positions of accreted subhaloes in the simulation, rather than use 
an analytic approximation. We do not attempt modifications to the 
models in this paper as our focus is on presenting and interpreting 
the observations. 

A fina l difference betw een this model and G14 is the adop¬ 
tion of the Marastonl 1 2005) SPS model, which differs from the 

y 


iBruzual A. & Charlotl ( 1993 ) model in its treatment of thermally- 
pulsating asymptotic giant branch stars (TP-AGB). These stars are 
important in the near-IR. The near-IR corresponds to the r and i- 
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Table 1. Summary of the different GALFORM models mentioned in this paper. 


Model 

Parent Model 

Cosmology 

Key Features 

Baueh et al. (2005) 

Cole et al. (2000) 

WMAPl 

Different IMFs in starbursting galaxies. 

Bower et al (2006) 

Cole et al (2000) 

WMAPl 

AGN feedback in hydrostatic haloes. 

Lagos et al (2012) 

Bower et al. (2006) 

WMAPl 

New star formation law accounting for 
fraction of hydrogen in moleculai* state. 

G14: Gonzalez-Perez et al (2014) 

Lagos (2012) 

WMAP7 

Updated cosmology. 

L14: Lacey et al (in prep.) 

Baugh et al (2005) & Lagos et al (2012) 

WMAP7 

Combines varying IMF 
and new star formation law. 


bands we use for luminosity and colour measurements iMarastonI 
l2005h . 

Note that the version of the L14 model we use is slightly dif¬ 
ferent to more recent versions of the Lacey et al (in prep.) model. 
The differences are tiny, however, and should not have any impact 
at the redshifts we consider. 

2.2.5 Model Stellar Masses 

We take stellar masses from the model, but rescale them hy h = 
0.702 (the model cosmology) in order to convert the physical units 
of the model, to the units of the observations, 

iMitchell et al.l (12013h demonstrates that applying broad-band SED 
fitting to model galaxy photometry can give a biased estimate of 
the true model galaxy’s mass. In order to investigate what effect 
this may have we show the inferred mass, from EqO against the 
true model mass in Eig.[T] In GI4 the relation between true and 
observed mass has a gradient, leading to smaller masses being sys¬ 
tematically under-predicted by the model. Red galaxies in G14 tend 
to be affected by this more than blue galaxies. In LI4 the differ¬ 
ences in mass have a larger scatter, the blue galaxies still show a 
gradient but the inferred mass is now generally an over-prediction 
of the true mass. The scatter between the inferred and true mass is 
also larger for the L14 model. Applying the colour relation to esti¬ 
mate mass, as compared to using the true mass, can cause both the 
masses and the typical colours of galaxies to change in the sample. 
ICampbell et alj 120151) demonstrate these differences can affect the 
clustering signal, with the effect being particularly marked for the 
L14 model. In this paper we take the model masses at face value, 
except for a small adjustment explained in Sec. 12.3.31 We discuss 
the possible effects of this when interpreting results. 

2.2.6 Literature GALFORM comparisons 

The GALFORM model has been compared to a variety of observa¬ 
tional measurements, including those made from this GAMA sam¬ 
ple. Indeed several of the GALFORM models were cre ated to ad¬ 
dress a disagreement with observa t ions. For e xample, iKim et alJ 
l2009h compared the ISower et"^ l2006l) and iFont et al.l l2008l) 
models to 2dFGRS data and found that the model predictions for 
the depende nce of clu s tering on luminosity were not successfully 
reproduced. iKim et al.l ( 120091) also found that the mod el had excess 
small scale clustering. The authors put this down to the lBower et alj 
( I 2 OO 6 I) model having to o many satellite galaxies. _ 

In support of this, iMcNaught-Roberts et alJ i20I4l) compares 
the luminos i ty fu nction as a function of environment of the 
iBower et alj ( l2006l) mocks to those of GAMA ga laxies, and finds 
an ex cess of faint, red galaxies in the model. iRobotham et alj 
( I 2 OI ih additionally found an excess of high-multiplicity (10 or 


( 5.8 



logio M*,true 


Figure 1. The inferred mass, from the Taylor et al (2011) colour to mass 
relation, versus the true mass of mock galaxies, for the G14 (top) and L14 
(bottom) models. The red and blue points represent red and blue galaxies, 
as selected by our colour cuts (Sec. 123). The contours are spaced evenly 
in log-space, between 5x10^ and 50x10^ galaxies dex~^. Infening the 
mass using obseivational relations, rather than using the hue model mass 
can make a significant difference and is in this context model dependent. 

more galaxies) groups in the [Bower et al.l ( l2006l) model compared 
to th e GAMA data. _ 

[Campbell et alj (1201. 5h combine SDSS and VIPERS measure¬ 
ments with the ones presented here to compare the clustering of 
galaxies as a function of stellar mass to that predicted by a variety 
of GALFORM versions, across several different epochs. They find 
good agreement between the models and data, but in order to fit 
small scale clustering they have to adopt a new, hybrid model for 
satellite orbits. In this model satellites follow resolved sub haloes 
and an analytic merger timescale is only computed once the sub¬ 
halo is lost. 

Their paper focuses on how different methods of estimating 
stellar mass affect clustering measurements; our paper takes the 
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complementary approach of comparing a wider range of galaxy 
pr operties. In addition , our pa per utilises the lightco ne modelling 
of iMerson et al.l ll201.3h . whilst l^mpbell et alj (1201.511 uses cluster¬ 
ing of galaxies in a single, constant cosmic time, snapshot. 

2.3 Sample Selection 

2.3.1 Redshift slices 

We want to study the evolution of galaxy clustering with luminos¬ 
ity, mass, colour and redshift. One approach to this is to use vol¬ 
ume limited samples, which are characterised by a uniform detec¬ 
tion probability across the sample volume. However, volume lim¬ 
ited samples reduce the amount of data available for the analyses, 
so we also consider magnitude-limited samples. In order to do this, 
we ensure the survey’s radial selection function is properly dealt 
with by the random catalogue (see Section IJTt . Using magnitude- 
limited samples means interpreting the results requires considera¬ 
tion of the selection function, but offers the most powerful test of 
the model. We separate all of our mass, luminosity and colour sam¬ 
ples into four redshift bins: low-z, 0.02 < z < 0.14, intermediate- 
z, 0.14 < 2 < 0.24, upper intermediate-z, 0.24 < 2 < 0.35 
and high-z, 0.35 < 2 < 0.5. The volumes of these slices are 

1.2 X 10® (Mpc//l)^ 5.0 X 10® (Mpc//l)^ 1.3 x lO'^ (Mpc/ft)® 
and 3.7 x 10^ (Mpc/fe)® respectively. 

In order to help differentiate between selection effects and 
galaxy formation effects, we also produce samples with upper red¬ 
shift cuts that have been decreased in order to make the sample es¬ 
sentially volume limited (comprise of at least 98% of galaxies with 
maximum observable redshifts greater than the sample limit )0 

2.3.2 Luminosity, mass and colour cuts 

Luminosity samples are produced from Eq. [T] with the measured 
SDSS r-band DR7 Petrosian magnitu de used to calculate t he ab¬ 
solute magnitude in the r-band, Af.. h . iLovedav et al.l (l2015h mea¬ 
sures the luminosity function, and finds the characteristic “knee” 
of the function at M* /, = — 20 . 6 ; past this magnitude the num¬ 
ber density of galaxies rapidly drops. This characteristic magnitude 
lies in our — 21.0 < Mr,h < — 20.0 sample, and we measure one 
sample of galaxies brighter than this and three fainter. 

Stellar mass is more complex to compute from observations, 
but is less complicated to predict from models. We use mass bins 
of 0.5 dex in size, except for the lowest mas s sample, where we 
increase the bin limits to have more galaxies. [Baldrv et al.l | |2012|) 
found for GAMA data at 2 < 0.06 that the characteristic knee of 
the mass function is at M* — lO'^® ®®M 0 /h^. This characteris¬ 
tic mass falls into our middle mass bin, we use two samples less 
massive and two samples more massive than this. 

In addition to luminosity and mass samples, galaxies are di¬ 
vided into a red population and a blue population. The colour of 
a galaxy is often used as a rough proxy for the age of its stellar 
population, with galaxies undergoing star formation generally be¬ 
ing expected to be bluer. It is therefore interesting to study how 

® Applying the same cuts to the mock catalogues leads to a larger frac¬ 
tion of some mock samples being magnitude limited, at worse the volume- 
limited fraction is 89%. Clustering from a mock catalogue for samples even 
more magnitude-limited than this (as low as 78%) show no significant dif¬ 
ferences with the clustering from the same mock with a deeper apparent 
magnitude cut (r < 21). We will therefore treat the mock measurements as 
volume limited. 


clustering differs as a function of colour. In Fig. we show the 
colour-magnitude diagram for our flux limited (r < 19.8) sample 
of galaxies. We use a sloping cu t with a gradient of 0.03, following 
the typical literature values (e.g. lSell et alj|2003l) . 

{g - r)o = -0.030 * {Mr,h - M*^h) + 0.678 (4) 

The intercept of this cut is such that there is an equal fraction of 
red and blue galaxies. We also show in Fig.|2]a volume-weighted 
colour-magnitude diagram. We see that the volume we observe has 
an obvious blue cloud, made up of faint galaxies, which is less ap¬ 
parent in the observed colour-magnitude diagram due to selection 
effects. The colour cut also separates the red sequence and blue 
cloud in the volume-weighted diagram. 

A summary of our cuts and the properties of our samples is 
given in Table 

2.3.3 Cuts and adjustments to the models 

When comparing the actual clustering of galaxies to models we 
want to uncover differences in how dark matter haloes are popu¬ 
lated. If we were to apply the same cuts on the models and the data 
we may get samples with very different number densities, since 
the luminosities or colours may be incorrectly predicted. Such dis¬ 
agreements are better explored by comparison of model predictions 
to other observational measurements like stellar mass functions, or 
colour distributions. In order to focus on the spatial distribution of 
comparable galaxies in the models and simulations, mock galaxy 
samples are often selected to have matching number densities to 
galaxies in the real comparison sample (e.e ISerlind et alj|2003l ; 
IZheng et al.ll2005| ; IContreras et alJl2013h . We therefore adjust the 
magnitudes of the mock galaxies in order for the luminosity func¬ 
tions of the mocks and data to match, in each of our four redshift 
slices. 

First we estimate the galaxy luminosity function of the r < 21 
unaltered mock magnitudes. We then reassign each galaxy in the 
mock a magnitude designed to map number densities in the un¬ 
altered mock luminosity function to the appropriate magnitude in 
the data luminosity function. We do this by computing the magni¬ 
tude difference between tbe cumulative mock luminosity function 
and the cumulative data luminosity function at the same cumulative 
number density. This forces the mock and data luminosity functions 
to agree. The faint end of the real luminosity functions suffer from 
incompleteness, as the different k-corrections of the galaxies re¬ 
sult in a colour dependence on which galaxies reach the apparent 
magnitude limit. In order to avoid this we estimate the complete lu¬ 
minosity function of real data in the three higher redshift slices by 
switching to the low-redshift real luminosity function past the point 
where colour-dependent completeness effects become important. 
The amplitude of the low-redshift luminosity function is scaled to 
ensure a smooth join with the higher redshift one. The mock cata¬ 
logue then has an r < 19.8 apparent magnitude limit applied to the 
shifted magnitudes. 

We apply the same k-correction polynomials to the mocks, but 
adjust the limits of the colour bins. This is to ensure the fraction 
of galaxies in each colour bin is the same as in the data. Thanks 
to this adjustment, applying the apparent magnitude limit results 
in a mock catalogue with close to the same selection function as 
the real data. This is important as many of our samples are magni¬ 
tude limited. Unfortunately, unless the mocks have the same colour- 
magnitude distribution as the data, it is impossible to perfectly 
model the selection function even with this technique of adjusting 
the mock k-corrections. This is because the magnitude distribution 

























Galaxy and Mass Assembly (GAMA): Projected Galaxy Clustering 1 


O 



O 


I 

02 






15-16-17-18-19-20-21-22-23 

tsj ui ^ o 

Ul Ul 




log 


lov 




15-16-17-18-19-20-21-22-23 

N) Ln '-J o 

L/1 1/1 


M°«-51ogio/i 


log 


lOV 


1.0 

0.8 

0.6 


o 





0.8 


0.6 


0.4 


0.0 


15-16-17-18-19-20-21-22-23 

0.0 c;i„„ j. ^Lnkjo 

ui Ln 

logio0 


M;-^-51ogio/i 




15-16-17-18-19-20-21-22-23 

NJ un '-J O 
L/i L/i 


-51ogio/i 


logio0 


0.6 


Figure 2. A restframe colour-magnitude diagram for the real GAMA galaxies and for the mock galaxies, as labelled. The top right panel gives the volume- 
weighted colour-magnitude diagram for the real GAMA galaxies. The red dashed lines shows our cuts to define red and blue samples of galaxies, this cut is 
different for each model and the data (see Section l2.3.3t . A histogram of colours has been added to the right of each panel. Contour levels are spaced linearly, 
the colour scale is the same in each of the colour-magnitude diagrams, except the volume weighted one. The contours were computed from a binned version 
of this plot, smoothed with a Gaussian filter. 


within a colour bin varies between real and mock data. The masses 
of the mock galaxies are rescaled to preserve the same mass-to- 
light ratio. Also note the mock catalogue is fit with the same P and 
Q as the data. 

We will now describe the cuts for colour selected samples. 
In Fig. [ 2 ] we can see that both models reproduce the bimodality 
of galaxy colours, but both predict a well-defined blue cloud not 
seen in the data. Note that whilst the data is flux limited, mean¬ 
ing faint galaxies may be lost, this selection effect is included in 
the model lightcone. One s uggestion of how to improve the mock 
galaxy colours was given in iFont et alj ( I 2 OO 8 I) : here the authors re¬ 


moved the unrealistic GALFORM assumption that all gas is removed 
from satellites as soon as they accrete onto a halo. By implement¬ 
ing a more realistic model with more gradual removal of gas by ram 
pressure stripping they found satellites could form stars for longer 
after being accreted, and ended up with colours between the red and 
blue sequence. However, the G14 and L14 models do not include 
any gradual ram pressure stripping. 

Looking at the models in more detail, we see the L14 has a red¬ 
der blue cloud than the G 14 model. This is expected as this model 
uses the iMarastorJ l l2005t) SPS model which has stronger near-IR 
emission for young stellar populations. We use the same gradient 
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as the data for the cut, and select the intercept in order to reproduce 
the observed fraction of red and blue galaxies. For the L14 model 
this intercept it 0.548, for the G14 model it is 0.498. From Fig.|2] 
we can see that the cut applied to the models clearly separates the 
red sequence from the blue-sequence mock galaxies. 

In Table IB^ and Table IB^ in Appendix lBl the sample sizes 
and properties for one realisation of GAMA are given for the two 
models. 

2.3.4 Comparison samples 

In addition to the samples we have describe d, we also produce a 
sample with which to compare our results to lZehavi et alj j201lh . 
In order to make the comparison samples as similar as possible, 
we use magnitud es corrected to = 0.1 and use the redshift 
cuts stipulated in IZehavi et alj ( l201lh . These magnit udes will be 
labelle d with the superscript ‘0.1’. Unfortunately, the lZehavi et alj 
( I 2 OI ih redshift cuts greatly restric t the v olume of our survey. As 
such, only two of the lzehavi et alj ( 1201 ih magnitude bin samples, 
-22.0 < < -21.0 and -23.0 < < -22.0 have a 

large enough volume in GAMA for worthwhile comparison, having 
volumes of 6.2x10® (Mpc//i)“® and 1.6x10® (Mpc//i)“® respec¬ 
tively. The redshift cuts for the fainter magnitude bins result in ex¬ 
tremely small volumes in GAMA (less than 1.1 x 10® (Mpc//i)“®). 

We also produce a series of volume-limited samples with 
which to compare to SDSS, using deeper redshift cuts appropri¬ 
ate to GAMA’s fainter selection limit. These samples are included 
in Table 

3 METHODOLOGY 

In this section we will introduce our method to produce random 
catalogues, before explaining how we compute the 2PCF, Wp(rp). 

3.1 Random catalogues 

3.1.1 The Cole (2011) method 

To measure clustering one needs a random set of points with the 
same radial and angular selection function as the data. We gen¬ 
et^ catalo gues of random positions using the method set out in 
ICold ( l201lh . which generates random catalogues from the real data 
in a wa y that remov es the effects of the large-scale structure. We 
use the Icold 1 I2OI ih method without the complications of the ex¬ 
tended method designed to fit a parameterised evolution model to 
the galaxies. For each galaxy in the catalogue the maximum vol¬ 
ume of space over which it could be observed, Umax, is calculated 
by finding 2min and 2max, the redshift where a galaxy meets the 
bright and faint magnitude limits of GAMA II. In addition to Umax, 
a density weighted maximum volume, Umax.dc, is calculated as 

rzmax ppp 

Umax.dc = / A{z)—dz (5) 

^min 

where A{z) is the over-density as a function of redshift and dV/dz 
is the comoving volume element per redshift element. Given these 
volumes, every real galaxy in the catalogue is cloned n times, 
where n is given by 

Unax 

72. = TlclonesT^ (O) 

►'max,dc 

with riciones being the total number of randoms divided by the num¬ 
ber of galaxies in the sample. Our default random catalogues have 


riciones = 400. When cloning the galaxy, all of the intrinsic galaxy 
properties are also cloned, so that the random points have a stel¬ 
lar mass, an absolute magnitude and a colour. The cloned galaxies 
are randomly distributed within the real galaxy’s Umax, with the 
GAMA angular mask used to ensure the angular selection function 
of the cloned galaxies matches that of the real galaxies. 

This method requires the estimation of A{z), which is done 
using an iterative method. Initially it is assumed A(z) = 1 every¬ 
where such that each galaxy is cloned the same number of times. 
From this random catalogue, A{z) is estimated from the redshift 
distribution of the randoms, nr{z), and the data, n^{z), using 

A{z) = riciones 

nr{z) 

A new random catalogue is then produced with this new estimate 
of A{z), and the whole process repeated until A{z) converges. We 
iterate 15 times, and find this is more than enough iterations for the 
process to converge. 

An added bonus of this method is tha t the Umax.d c can be used 
to easily estimate the luminosity function ( IColell20flh . We use this 
technique when estimating luminosity functions for the purpose of 
adjusting mock galaxy magnitudes to match the abundance of real 
galaxies. A study of the lumi nosity functions of G AMA galaxies 
from this estimator is given in lLovedav et al.l ( l2015h . 

3.1.2 Windowed Clones 

As previously mentioned, the creation of this random catalogue re¬ 
quires that one adopt a certain value of Q (Eq.[T]l. Additionally, our 
choice of Q is motivated by our decision not to consider evolution 
in the overall number density of galaxies. It may be the case that 
our chosen Q is not true to the data, or some density evolution is 
present in the sample. To mi t igate p ossible effects from this we add 
a new technique to the lColel 220111) method. 

The key idea of this method is to restrict the redshift of cloned 
galaxies to some window function around the redshift of the origi¬ 
nal galaxy. Evolution is therefore included naturally in the random 
catalogue, as cloned galaxies are kept close to their original red¬ 
shift. Hence the error introduced by the inadequacy of the adopted 
P and Q model is limited. We choose to define our window as a 
function of volume, W(V); this is because volume is most closely 
related to the expected fluctuations in the galaxy A{z). With a win¬ 
dow defined in terms of volume, the n{z) is smeared out more in 
the sample variance prone, low-redshift part of the n{z) and less 
at higher redshift. We adopt a Gaussian window function, with 
a = 3.5 X 10®(Mpc/Ii)®, truncated at 2a. We found this to be 
a good compromise between limiting the effects of evolution and 
smoothing out large-scale structure. We also include this window in 
the computation of Umax and Umax.dc, by weighting volume slices 
by the window when numerically computing the integral. When the 
window function reaches a boundary, either the Umax of an object 
or the limit of the survey, it is reflected from that boundary in or¬ 
der to ensure the randoms are scattered symmetrically across the 
volume. Simply truncating the window at a redshift limit would 
result in randoms being artificially more likely to be moved away 
from the limit. This method will be further explored, along with 
the choices of a, in a future paper. We will, however, demonstrate 
one success of the method here empirically. In Eig.j^we show the 
redshift distribution of GAMA cloned galaxies, for different bins 
of the original galaxy redshift. In the top panel of Eig.[3]we can see 
that cloned galaxies are spread widely across the full redshift range 
of the survey. In the bottom panel of Eig. [3 we can see that the 
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Table 2. Different galaxy samples, sample sizes and median properties. The subscript ‘med’ indicates the values are medians. Samples with redshift limits 
marked with asterisks are magnitude limited, so should be treated with careful consideration of the GAMA selection function. In samples without an asterisk 
at least 98% of the members are volume limited. Values in brackets are rms scatter. 


Sample 

'^min 

-2^max 

-^gals 

•^med 

AT’med 

loglo(^*/M0/l ^)med 

id ^)o,med 

-18.00 < Mr,h < -17.00 





-17.46 (0.28) 

8.62 (0.29) 

0.42 (0.20) 



0.14* 



-17.62 (0.27) 

8.67 (0.33) 

0.41 (0.19) 

-19.00 < Mr,h < -18.00 

0.02 

0.11 

6950 

0.09 

-18.47 (0.29) 

9.13 (0.31) 

0.47 (0.21) 


0.02 

0.14* 

13149 

0.11 

-18.54(0.27) 

9.15 (0.31) 

0.47(0.19) 


0.14 

0.24* 

3307 

0.15 

-18.85(0.13) 

9.26 (0.28) 

0.46 (0.16) 

-20.00 < Mr,h < -19.00 

0.02 


11741 

0.11 

-19.47(0.29) 

9.74 (0.29) 

0.62(0.17) 


0.14 


7801 

0.16 

-19.46 (0.29) 

9.74 (0.31) 

0.60 (0.22) 


0.14 


27755 

0.18 

-19.59(0.27) 

9.79 (0.32) 

0.59 (0.20) 


0.24 


2834 

0.26 

-19.90(0.10) 

9.81 (0.38) 

0.45 (0.18) 

-21.00 < Mr,h < -20.00 


0.14 

6945 

0.11 

-20.40 (0.28) 

10.24 (0.25) 

0.71 (0.15) 


0.14 

0.24 

22327 

0.20 

-20.39 (0.28) 

10.27 (0.28) 

0.70 (0.15) 


0.24 

0.35* 

37541 

0.28 

-20.55 (0.27) 

10.39 (0.28) 

0.68(0.17) 


0.35 

0.50* 

3593 

0.37 

-20.86 (0.14) 

10.42 (0.28) 

0.55 (0.16) 

-22.00 < Mr,h < -21.00 


0.14 

1843 

0.12 

-21.27(0.25) 

10.67 (0.19) 

0.74(0.11) 


0.14 

0.24 

5865 


-21.28 (0.25) 

10.70 (0.21) 

0.74 (0.13) 


0.24 

0.35 

15353 

^^9 

-21.29 (0.25) 

10.78 (0.23) 

0.74 (0.17) 


0.35 

0.37 

3168 

0.36 

-21.28 (0.25) 

10.82 (0.23) 

0.74(0.14) 


0.35 


16114 


-21.42 (0.27) 

10.86 (0.28) 

0.73 (0.22) 

-23.00 < i < -22.00 

0.01 

0.50 

2273 

0.37 

-22.25(0.19) 

11.26 (0.28) 

0.79 (0.18) 

-22.00 < Mjj < -21.00 

0.01 

0.38 

23746 

0.29 

-21.36 (0.25) 

10.79 (0.23) 

0.74(0.15) 

-21.00 < < -20.00 

0.01 

0.26 

33965 

0.20 

-20.47 (0.28) 

10.31 (0.28) 

0.71 (0.16) 

-20.00 < Mjj < -19.00 

0.01 

0.18 

22487 

0.14 

-19.53(0.29) 

9.78 (0.31) 

0.62(0.19) 

-19.00 < < -18.00 

0.01 

0.12 

8129 

0.09 

-18.52(0.29) 

9.16(0.31) 

0.48 (0.21) 

8.50 < logjQ Mt,/MQh~^ < 9.5 

0.02 

0.05 

1049 


-17.95 (0.72) 

8.88 (0.28) 

0.46 (0.15) 


0.02 

0.14* 

18533 


-18.42(0.59) 

9.07 (0.27) 

0.44(0.13) 


0.14 

0.24* 

8472 


-19.19(0.33) 

9.32 (0.16) 

0.39(0.10) 

9.50 < logio M*/M0 /i-2 < 10.00 

0.02 

0.14 

10053 

0.11 

-19.43(0.47) 

9.74 (0.14) 

0.66 (0.12) 


0.14 

0.24* 

19173 

0.18 

-19.64(0.37) 

9.79 (0.14) 

0.56(0.12) 


0.24 

0.35* 

5331 

0.27 

-20.11 (0.25) 

9.84 (0.13) 

0.41 (0.11) 

10.00 < logio M*/M0h-2 < 10.50 

0.02 


7385 

0.11 

-20.28 (0.44) 

10.21 (0.14) 

0.73 (0.12) 


0.14 

0.18 

7749 

0.16 

-20.25 (0.43) 

10.21 (0.14) 

0.73 (0.09) 


0.14 

0.24* 

23216 

0.20 

-20.24 (0.40) 

10.22 (0.14) 

0.72(0.10) 


0.24 

0.35* 

23584 

0.28 

-20.48 (0.29) 

10.32 (0.14) 

0.63(0.12) 


0.35 

0.50* 

3132 

0.38 

-20.92 (0.26) 

10.36 (0.13) 

0.48(0.10) 

10.50 < logio M*/M0h-2 < 11.00 

0.02 

0.14 

2189 

0.12 

-21.15(0.38) 

10.65 (0.12) 

0.76(0.10) 


0.14 

0.24 

8018 

0.20 

-21.10(0.38) 

10.65 (0.13) 

0.76 (0.09) 


0.24 

0.29 

9968 

0.27 

-21.00 (0.38) 

10.66 (0.13) 

0.76 (0.09) 


0.24 

0.35* 

24416 

0.30 

-21.02 (0.35) 

10.68 (0.13) 

0.75 (0.09) 


0.35 

0.50* 

11800 

0.38 

-21.31 (0.26) 

10.78 (0.14) 

0.69 (0.12) 

11.00 < logio M*/M0h-2 < 11.50 

0.24 

0.35 

2572 

0.31 

-21.86(0.38) 

11.10(0.10) 

0.80(0.15) 



0.37 

698 

0.36 

-21.74 (0.35) 

11.11 (0.11) 

0.82 (0.12) 



0.50* 

5962 

0.41 

-21.81 (0.29) 

11.14(0.12) 

0.82(0.11) 

Red 

0.02 

0.14* 

15842 

0.11 

-19.61 (1.12) 

9.98 (0.54) 

0.73 (0.17) 

(g - r)o + 0.03iMr,h - M*^) > 0.678 

0.14 

0.24* 

29686 

0.19 

-20.16(0.70) 

10.26 (0.35) 

0.75 (0.14) 


0.24 

0.35* 

30465 

0.29 

-20.81 (0.50) 

10.62 (0.29) 

0.76 (0.15) 


0.35 

0.50* 

11461 

0.39 

-21.51 (0.39) 

11.00(0.29) 

0.80 (0.23) 

Blue 


0.14* 

26437 

0.10 

-18.61 (1.21) 

9.11 (0.62) 

0.44(0.12) 

(g - r)o + 0m{Mr,h - M*^) < 0.678 


0.24* 

29906 

0.19 

-19.77 (0.67) 

9.74 (0.42) 

0.51 (0.11) 



0.35* 

26227 

0.29 

-20.56 (0.49) 

10.26 (0.35) 

0.56 (0.12) 



0.50* 

10122 

0.39 

-21.22 (0.39) 

10.61 (0.30) 

0.58 (0.12) 
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Figure 3. The distribution of the final redshifts of the random points, in 
bins of their initial redshift, for a random catalogue for the real GAMA data, 
generated without (top) and with (bottom) the window function update. The 
window function update acts to limit how far a cloned galaxy can move from 
its initial position. This in turn limits the effects of any unmodelled galaxy 
evolution on the random catalogue. 

addition of a window function has limited the redshift range over 
which a cloned galaxy can move. As such it has limited the effects 
of deviations of the evolution away from the description given by 
the adopted values of P and Q. 

This windowing method requires further testing and tuning, to 
see how successful it is in removing the effects on the random cata¬ 
logue from unmodelled galaxy evolution. However, for this work 
we show in Appendix that the windowing affects our results 
much less than sample variance. 

3.1.3 The resultant catalogue 

We show in Fig. |4] (top) the redshift distribution of the data and 
the randoms for different iterations of this process. We see that us¬ 
ing the density-corrected maximum volume (olive and blue dashed 


line) only introduces subtle differences into the ni{z) of the ran¬ 
doms, as compared to simply using Vmax (red dashed line). The 
twelfth iteration (blue) and the fifteenth iteration (olive) agree; this 
indicates the process has converged. The nr{z) of the randoms is 
a good fit to the ng(z) of the total sample; later we will check 
the random nr(z) is appropriate for galaxies split into magnitude, 
colour and mass samples. Note that it does appear that the random 
catalogue slightly follows the under-density at ^ ~ 0.22; we found 
this was not the case when not using the window function and as 
such claim this is an unfo rtunate side effect of the modification to 
the method (lFarrowl2013l) . We have, however, tested our clustering 
results using randoms from the old and new methods, and only find 
small differences that are at large scales. These differences are a 
fraction of the error bars, and as such our results will not be signif¬ 
icantly affected. 

In Fig.|4](Bottom) we plot the over-density estimate, Eq.|71 for 
successive iterations of our random catalogue generating method. 
As expected, the iterations ac t to slightly increase the over-density 
estimates, as the ICold (1201 ih method acts to remove their effect 
from the random catalogue. This method of generating randoms 
ameliorates the effects of large-scale structure on the generation of 
a random catalogue, as over-represented galaxies from over-dense 
regions are cloned fe wer times wh ilst under-represented galaxies 
are cloned more times. IColel(l201ll) demonstrates how this scheme 
can produce a random catalogue unbiased by large-scale structure. 
A strength of this approach to generating randoms is that the ran¬ 
dom catalogue comes with all of the properties of the galaxy cata¬ 
logue. One can then apply the same selection to the random cata¬ 
logues and the galaxy catalogues so that the random catalogue has 
the correct angular and radial distribution. We intend to make these 
random catalogues available to the GAMA community via a DMU, 
which we also intend to distribute to the larger astronomical com¬ 
munity once the GAMA data become public. 

In Fig.[5]we show the redshift distribution of the randoms and 
GAMA data, split into samples. We see for the luminosity, mass 
and colour samples we study the random redshift distribution is an 
excellent fit to the data. 


3.1.4 Randoms for the Mocks 

We apply the same techniques to generate randoms for the mock 
catalogues. This technique is applied first to the original mock cat¬ 
alogue with unadjusted magnitudes and r < 21.0, in order to esti¬ 
mate the mock luminosity function. The technique is applied again 
to the mock catalogue with adjusted magnitudes and r < 19.8, to 
yield the final random catalogue. We find that the resultant n{z) of 
the randoms are a good fit to the mock catalogues. 


3.2 Projected clustering 

We measure t he 2PCF of galax ies using pairs of galaxies and ran¬ 
doms with the iHamiltonl d 19931) estimat or. Following th e standard 
approach adopted in the literature (e.g. ICoil et al]l2008l) we mea¬ 
sure pair separations parallel, tt, and transverse, Vp, to the line of 
sight for each pair. These are computed by first converting the an¬ 
gular position and redshift of each object to a vector, r. We then 
define a line-of-sight direction to a pair as 1 = (ri + r^)/2, where 
ri and r 2 are the positions of the two pair members. The parallel 
to the line-of-sight distance, tt, is the projection of the separation. 
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Figure 4. (Top) The redshift distribution of GAMA II galaxies (solid) and our randoms (dashed) for multiple iterations of the Cole (2011) random catalogue 
generating approach (see legend), combined with our window function modifications, as explained in Section 4.3.2. The redshift distribution of the randoms 
is a good match to the data. (Bottom) Our estimates of the galaxy over-density as a function of redshift, from the ratio of the galaxy and random redshift 
distributions. The error bars are from the rms scatter of the 26 mock realisations, as such they account for the large-scale structure variations. Note these error 
bars may be underestimates, as the mock realisations over-sample the simulation volume at high redshifts. 


s = r 2 — ri, onto the line of sight 

s - 1 



( 8 ) 


The separation transverse to the line of sight is then 

Tp = ^|sP - TT^. (9) 

Pairs are binned onto a grid of tt and Vp. In order to maximise 
signal-to-noise ratio we use a mixed linear and logarithmic bin¬ 
ning scheme. Bins are linear with 0.0625 /i“^Mpc width up to 
0.3 /i“^Mpc, and are 0.12 dex in size for bins larger than that. 
Whilst we generate 400 times as many random points as the data, 
when computing clustering we use between Udones = 32 and 
ticiones = 200 times more randoms than data, the number varies in 
order to maintain good statistics for the samples where the galaxies 
only sparsely populate the volume. 

It is well known that the measured correlation function is dis¬ 
torted by the peculiar velocities of galaxies. On larger scales the 
infall of galaxies squ ash the obser ved correlation function in the 
line-of-sight direction iKaiseJl987h . On small scales the virial mo¬ 
tions of galaxies within clu sters can elong ate the correlation func¬ 
tion along the line of sight l ljacksotj|l972h . These distortions have 
been studied in th e GAMA data by Loveday et al (in prep.) and 
[Blake etalJj2013l) ; in this work we focus instead on the projected 
correlation function, Wp(rp). The projected correlation function 


is a standard approach for dealing with redshift space distortions, 
which involves integrating ^(rp, tt) along the tt direction to min¬ 
imise their effects, thus 

Wp(rp) = 2 / ^(rp,7r)d7r. (10) 

Jo 

In practice this integral is carried out numerically using our 
^(rp, 7 r) grid. The choice of TTmax warrants careful consideration. 
Ideally, one would use the largest possible value of TTmax to in¬ 
clude most of the 2D clustering signal, and because theoretically 
the effects of redshift space distortions are only removed if you in¬ 
tegrate out to infinity. Unfortunately, in real surveys noise affects 
the measurements, and measurements for large values of TTmax can 
be particularly noisy. We adopt a value of TTmax = 41Mpc/i“^ 
for our measurements, which is a reasonabl e TTmax value to con - 
sider fo r the scales analysed (see e.g. Fig. 1 of iNorberg et al.l2009h . 
For the IZehavi et alj ( 1201 ih comparison sample we use TTmax = 
60 Mpc h~^, which m akes our measurements noisier but matches 
the lZehavi et alJ ( 1201 ih TTmax- 

The results of Eq.[T0l with TTmax = oo, can be calculated an¬ 
alytically for a spherically symmetric power law correlation func¬ 
tion, ^(rp) = (rp/ro)~^, where ro and 7 are constants. The result 
is 


Wp(rp) =rp ( — 
. '1 


roVr(l/2)r((7-l)/2) 


r(7/2) 


( 11 ) 
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where F is the Gamma Function. We will fit Eq.[TT]to some of our 
samples in Section [5] in order to measure correlation lengths, ro. 
As a power law is not a good fit over the whole 2PCF we restrict 
the fit to the scales 0.2Mpc/i~^ < rp < 9.0 MpcTo fit 
we adopt a least-squares minimisation method using the diagonal 
terms of the covariance matrix. To check how the variation of 7 
could influence the derived value of ro, we also fit power laws with 
a fixed 7 = —1.8. All of the measured values are given in the 
Appendix in Table I Fill We see only a small ditference between the 
best fitting ro values where 7 is free to vary and where 7 = — 1 . 8 . 
The difference is not large enough to affect our conclusions. 

When plotting clustering, we often include a reference power 
law line or divide through by this reference power law, Wref, to 
allow easier comparison between plots. We use the IZehavi et alj 
(l 201 lh power law fit to their — 21.0 < < — 20.0 sample for 

this purpose, which has ro = 5.33/i“^ Mpc and 7 = —1.81. 

The clustering results in this paper have been cross-checked 
numerous times to great accuracy against independent clustering 
analyses. 


Ie3 



le3 



Figure 5. The redshift distribution of the data (solid lines) compared to the 
redshift distribution of the randoms (dashed lines), for different samples 
as indicated in the legend. Vertical dotted lines mark the positions of our 
redshift cuts. The randoms provide an excellent description of the data. 


3.3 Error estimates 

To compute error bars on our cl ustering measureme nts for GAMA 
data, 9 jack-knife samples (e.g. IZehavi et alj| 2002 h . 3 per region, 
are formed by rejecting roughly equal-area regions of data. From 
this method, the covariance matrix Cij, is calculated, the square- 
ro ot of the diagona l terms of which give the error bars. Work such 
as iGuo et alJOOldl) show that with large surveys like SDSS, a large 
number of jack-knife resamplings (~ 100 ) can give reliable esti¬ 
mates of the covariance matrix. Our smaller area, necessitating a 
smaller number of jack-knife regions, means we should further test 
our covariance matrix estimates. 

We also test how unbiased using jack-knife errors is for our 
sample. For this we utilise the 26 realisations of each model, cre¬ 
ated by considering different lightcones through the simulation. 
These estimates of the error should be more realistic, as they are 
built from a series of lightcones which together sample a larger vol¬ 
ume than the GAMA data, even though they are drawn from a sin¬ 
gle simulation. This latter fact implies that the mock errors shown 
are likely underestimates of the true sample variance, as estimated 
from an series of lightcones constructed from independent simu¬ 
lations. We computed the jack-knife errors, for a single mock, and 
the root mean square (rms) scatter between the different realisations 
for all of the samples selected on mass. In Fig. we show the ra¬ 
tio of the two error estimates as a function of projected separation. 
Different colours are from the different models, as indicated. The 
dashed and solid lines represent two different realisations of the 
model, encompassing different volumes of the simulation. These 
two volumes are the same for both models, in order to disentangle 
the effects of sample variance and galaxy formation. 

In the highest redshift slices the error estimates converge, 
however note that this convergence could be due to the two realisa¬ 
tions sampling some of the same volume of the N-body simulation 
(see Section l2.2.2t . In the two lower redshift slices the realisations 
should be nearly independent. We can see that for some redshift and 
mass ranges the jack-knife errors give the same uncertainty as the 
model, whilst for others the errors disagree. In particular, the low- 
redshift jack-knife errors give a lower estimate in the first volume. 
This suggests for our samples jack-knife errors could be underes¬ 
timates at low redshift. As expected, we only see a small differ¬ 
ence between the ratios for the two models in the same volume, as 
such it appears that the ratio is most sensitive to sample variance. 
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Figure 6. The ratio of error computed from the scatter between 26 mock realisations, and jack-knife error estimates from two different realisations. Red lines 
are from the L14 model whilst blue lines are from the G14 model. We can see the en'or estimates from jack-knife resampling can be very different to those 
from the mock realisations. 


Despite the limitations of the jack-knife error estimates with our 
samples, we use them with the data as it has the advantage of not 
being model dependent, nor sensitive to how accurately the models 
reproduce the observed GAMA clustering. However, for the mock 
predictions we will always plot the error computed from the scatter 
between the mocks. This is not only true for the clustering mea¬ 
surement, but also for the errors on the parameters of our powerlaw 
fits. When we plot the mean correlation functions of our 26 mock 
realisations, we do not divide our error estimates by y/^, such that 
the errors reflect the scatter expected on an individual realisation of 
the GAMA survey. Again recall these errors are underestimates in 
the two higher redshift slices due to the overlap between lightcones 
of different realisations in the simulation. 


When assessing the goodness-of-fit of the mock predictions to 
the data, or our data to literature data, we will use the covariance 
matrix. For the comparisonsof our GAMA observations to the liter¬ 
ature we use the covariance matrix from the ja ck-knife resampling . 
As we compare to the SDSS measurements of IZehavi et al.l | |201 ih . 
we can disregard the contribution of the errors on the literature mea¬ 
surements to the goodness-of-fit as they come from a much larger 
volume than our GAMA samples. When comparing the mocks to 
the data, we use the covariance matrix computed from the 26 mock 
realisations. This tests the hypothesis that the data is a realisation of 
the mock. As with the SDSS measurements, the error on the mean 


measurements from the 26 mocks are small compared the real data 
uncertainties and are therefore disregarded. 

The computation of the inverse of the covariance matrix, 
can cause problems as our estimates of the covariance, from 
a limited number of resamplings or realisations, can be noisy. 
Our m etho d to correct for this follows iGaztanaga & Scoccimarrd 
(l2005h and iMarm et alj ( l2013ll . We first compute the correlation 
matrix. 


Ci, 




( 12 ) 


where (Ji is the standard deviation of the 2PCF measurement 
and i and j are indices running over all of the 2PCF measure¬ 
ments. We then carry out a singular value decomposition (SVD) 
of it, yielding C = UCsvdU^, where U and are rotation 
matrices and Csvd is a diagonal matrix with elements XijSij. The 
rotation matrix acts to transform the data points into a coordinate 
system where they are no longer correlated. The basis of this new 
coordinate system are the eigenvectors of the covariance matrix, 
given by the columns of the matrix U, which have XiAg as t heir 
eigenvalues. As explained in iGaztafiaga & Scoccimarrd ( l2005h the 
eigenmodes of the data expressed in this basis should be uncorre¬ 
lated, Gaussian distributed and given by 


W(i} = 


WU) 


(13) 
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where W(j) are elements of a vector of the 2PCF measurements, 
and W (i) are the new eigenmodes. As we have a noisy estimate 
of the covariance matrix, from a limited number of jack-knife sam¬ 
ples or mock realisations, some of the eigenvectors will be poorly 
estimated and using them would bias the estimated x^- We would 
expect the modes which contribute least to the the variance to be 
most likely to suffer from this problem; these modes have small 
eigenvalues. 

We _ experiment ed _ with applying the 

iGaztanaga & Scoccimarrol (l2005h cut of \ij > ^2/Wmock,jk 
where Wmockjk is either the number of mock realisations for the 
tests of the model (Wmockjk = 26), or the number of jack-knife 
resamplings for th e x^ tests of whether our results agree with 
IZehavi et al.l (l201lh (Wmockjk = 9). However, we found that this 
could, counter intuitively, lead to noisier data sets or data sets with 
fewer reliable measurements having more accepted modes. We 
therefore adopted the approach of only taking the 4 largest modes, 
which for our data is sligh t ly mo re conservative than applying the 
IGaztanaga & Scoccimarrol ^20051) cuts. 

To invert Csvd we simply take 1 / Xij for each element, this is 
correct as the matrix is diagonal. We set 1 /Xij — 0 for eigenvalues 
failing the cuts. Setting values to zero in this way means some pos¬ 
sible degrees of freedom are removed. We compute x^ using the 
deviations divided by their associated errors, i.e. 

2 ^ (H/(i)data - fy Wref) (W"0')data " VKO^ef) 

X — 

(Ti Ci 

(14) 

Here the ‘data’ measurements are the measured GAMA 2PCF and 
the ‘ref’ measurements are either literature data or the 2PCF predic¬ 
tion from the combined mock catalogues. As our x^ values are still 
likely to be somewhat inaccurate, we will only indicate whether the 
cumulative x^ distribution for four degrees of freedom suggests a 
probability of a certain measurement is less than 2 %, demonstrat¬ 
ing it shows highly statistically significant differences. 


4 RESULTS 

4.1 Comparison to literature results 

In this section w e compare our measurements to those of 
IZehavi et al.l (1201 ih . which come from the SDSS. This comparison 
acts as a test that our methods for fc-correcting galaxies, producing 
randoms and measuring clustering give reasonable results. We can 
also use the large area of the SDSS to gauge if the GAMA volume 
is particularly under-dense or over-dense. 

In Fig. 7] we s how our measurements (blue) against those of 
IZehavi et al. OOl ih (black), both for the same redshift cuts. We also 
show, in red, our measurements for deeper volume limited samples 
using the GAMA apparent magnitude limit in Fig.|7] 

When the p robability of the hy pothesis that these measure¬ 
ments agree with IZehavi et alj (1201 11 ) falls below 2% i t is indicated 
on the panel. As a simplification we do not use the IZehavi et alJ 
( I 2 OIII) error estimates in this calculation, as our errors are the dom¬ 
inant source of uncertainty. The samples with the SDSS redshift 
limits show good agreement with the SDSS measurements, sug¬ 
gesting the GAMA volume is not particularly unusual in the red¬ 
shift ranges probed by these samples. 

The major ity of our deeper, G AMA volume-limited samples 
agree with the IZehavi et alJ ( I 2 OI ih measurements. The brightest 
how ever shows a s ignific ant variation compared to the lower red¬ 
shift IZehavi et al.l ( 1201 ih measurements. Note, however, that the 


jack-knife covariance matrices may not be accurate for a sample 
with such a small size and large Poisson errors. The faintest sam¬ 
ple also appears to have a highly st atistically significant difference 
compared to the lZehavi et alJdTOllh measurement. Specifically, the 
clustering signal has a lower amplitude. Note, however, the vol¬ 
ume of this sample is even smaller than our low -2 slice, for which 
Fig.l^demonstrates the jack-knife errors can be underestimates by 
a factor of ^ 3. In support of this being just being due to sample 
variance [Driver et ^ ( 1201 ih find that the GAMA survey is 15% 
under-dense compared to SDSS DR7 up to 2 = 0.1, and the me¬ 
dian redshift of this sample is 2 = 0.09. 


4.2 Stellar mass dependent clustering 

4.2.1 The full shape of the correlation functions 

We will now begin to study clustering as a function of a quantity 
that requires less modelling in the mocks but more modelling in the 
data: stellar mass. In Fig. we show the clustering of galaxies as a 
function of mass, divided by our reference power law. Dividing by a 
reference power law in this way means our measurements show the 
square of the galaxy bias, relative to the fiducial power law. Note 
that not all of the samples are volume limited, so one should be 
careful to use Tablel^to characterise the typical masses, luminosi¬ 
ties and colours of the samples. Samples that are volume limited 
are indicated by a star in Table|2 Samples which are more effected 
by the magnitude selecti on are expected to hav e less clustering than 
volume-limited samples dMeneux et al.ll20()8h . and are expected to 
be bluer due to our colour-dependent fe-corrections. This can be 
seen in Table |2 where the higher redshift, more magnitude limited, 
samples of a particular magnitude or mass range tend to have bluer 
median colours. The mocks have been constructed to have the same 
selection function, so it is fair to compare directly. 

In all the redshift slices we observe more clustering in the 
more massive galaxy samples. This is only seen in Fig. [ 8 ] with 
volume-limited samples in the lowest redshift slice. However, in 
the two highest redshift slices the median colours and magnitudes 
of the second most massive sample, which is magnitude limited, 
are fairly similar to those of the low redshift volume-limited sample 
(Table O. This similarity suggests that the large change in cluster¬ 
ing amplitude seen between the 10 . 5 < logj^QM*/M 0 fe ^ < 11.0 
sample and the 11.0 < log^Q M*/M 0 /i“^ < 11.5 sample is real, 
even though the less massive of these samples is magnitude-limited. 
In the redshift slices we also see the small-scale clustering (around 
< 3 /i“^Mpc), relative to the large-scale clustering, increases with 
mass. This could be interpreted as an increasing fraction of satel¬ 
lite galaxies in the higher mass samples, or it could be a result of 
fainter, redder satellite galaxies being lost in the more magnitude- 
limited samples. 

The mock predictions for the two models are fairl y similar to 
one another. This similarity suggests the adoption of the ijiang et all 
(|2008L|2 oI^ approach to computing satellite merger times has only 
a small effect on the clustering. The general trend of more massive 
galaxies being more clustered is reproduced. In a ACDM cosmol¬ 
ogy this implies in both the model and the data more massive galax¬ 
ies reside in more massive dark matter haloes. As discussed, some 
of this effect may also be due to using magnitude-limited rather 
than volume-limited samples. While the amplitude of the cluster¬ 
ing in the models is generally an acceptable match to the data, the 
small-scale clustering is often incorrectly predicted. 

For the two slices below 2 = 0.24, galaxy samples in the 
10.0 < logj^Q M*/M 0 fe“^ < 10.5 mass range appear, by eye, to 
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Figure 7. Clustering measurements from Zehavi et al. (2011) (black), along with our measurement of clustering for the same magnitude and redshift cuts (blue) 
and for our, deeper volume-limited samples (red). Error bars are from jack-knife resamplings of the data. As GAMA has a smaller area than SDSS, applying 
the same magnitude and redshift cuts as Zehavi et al. (2011) resulted in the fainter samples being too small to measure clustering. Where the probability of 
our data being a realisation of the Zehavi et al. (2011) data falls to less than 2%, this is indicated on the panel with the label P(x^) < 2%, colour coded by 
the sample tested. The measurements using the same redshift limits as Zehavi et al (2011) all pass this criteria. 


have too high small-scale clustering. This discrepancy is borne out 
statistically, with values less than 2 % probability of this data being 
a realisation of the model. For the G14 model this over-prediction 
of small scale clustering persists to higher redshift. The mock pre¬ 
dictions for the 10.0 < logj^g < 10.5 sample at 

0.35 < z < 0.5 also show a significant discrepancy, which does 
not appear to be restricted to small scales. However recall that the 
covariance matrices derived from the mocks are underestimates of 
the true error in the high redshift slice, due to the oversampling of 
the N-body simulation. 

Additionally, the 8.5 < logj^g < 9.5 sample in 

the 0.14 < 2 < 0.24 slice and the 9.5 < logj^g < 10 

in the 0.24 < 2 < 0.35 slice has a clustering signal on scales 
less than a few Mpcs that is significantly too small compared to 
the GAMA data. This is true even for the latter of those samples, 
where the data error bars might suggest the model is a good fit. This 
is because the error bars are much smaller in the mock predictions 
than in the real data, as the mock number densities are much larger 
for this sample (recall only the luminosity functions were forced to 
match). 

In Fig. 13 we plot samples which have had their upper redshift 
limits adjusted in order to make them volume limited. In general 
they imply the same conclusions as we drew from the flux-limited 


samples. We again see that the more massive galaxy samples are 
more clustered. Several of the mock samples have too much clus¬ 
tering on small scales. Seeing this effect in the volume-limited sam¬ 
ples means it is not related to the modeling of the selecton function. 

4.2.2 Results from power law fits 

To further investigate the redshift and mass evolution of the clus¬ 
tering, we show in Fig. [T3 the rg values from the power-law fits 
to the measurements, as a function of the ratio of the median sam¬ 
ple mass to the characteri stic mass of the ste llar mass function at 
2 = 0.0, lO^°'®®M 0 /h^ tealdrv et al.ll2012h . Samples consisting 
of 95% or more volume-limited galaxies are indicated by a star, 
magnitude limited samples are indicated by a triangle. We see very 
similar dependence of clustering on mass in both the GAMA data 
and the models. Also, in the magnitude-limited samples of the data 
and the models the increase in clustering strength as a function of 
stellar mass is faster for the higher redshift samples. This effect is 
also seen in the mock catalogues. This can also be viewed as less 
massive samples evolving faster with redshift than more massive 
ones. However, note this effect is not seen when considering the 
volume-limited samples, which only show weak evidence of red¬ 
shift evolution. 
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Figure 8. The projected two-point correlation function of real GAMA (black), G14 model (red) and L14 model (blue) galaxies as a function of redshift 
(different columns) and stellar mass (different rows). The measurements are divided by the reference power law defined in Section [Til EiTor bars are from 
jack-knife resamplings for the real GAMA data and the scatter between the mock realisations for the mocks. Marked with P < 2% are samples for which the 
hypothesis that the data is a realisation of the model has less than a 2% probability, colour coded according to the model. Despite clustering measurements 
not being used when creating the models, they reproduce the same trends with mass and redshift as the data. Some samples, however, do show significant 
differences, particulai'ly on small scales. 


In order to further explore possible red shift evo l ution, in 
Fig.|T0|we plot the bias versus mass relation oflEi et alJ (l2006h as 
dashed lines, colour-coded according to redshift. We convert this 
relation into ro values assuming a power-law correlation function, 
and calculating the ro of dark matter using the bias and ro of the 
galaxy sample containing M*, i.e. ro_dm = ro(M*)6(M*)“^/^. 
We can then use the relation ro{M) — (b(M))y^ro.dm . We 
use our own power law fit to the published iLi et al.l ( l2006ll sam¬ 
ple containing M* for ro(M*l . A co rresponding bias, b{M*), is 
taken directly from the lLi et Tl 1 I 2 OO 6 I) fitting formula. The dashed 
lines at differ ent redshift s are calculated using the passive evolu¬ 
tion model of iFrvI il99(i) to evolve the bias values, and using the 
growth factor D{z) to evolve the dark matter ro as ro,dm(z) = 
''^o,din{z'){D{z)/D{z'))'^^~'. Such a model gives the expected evo¬ 
lution of ro for a model where galaxies formed in some density field 
before moving along trajectories at a velocity defined by their lo¬ 
cal gravitational potential. The growth factor was calculated for ou r 
cosmology from the approximate formula of ICarroll et alj jl992h . 
We assume 7 = 1.8 for the power laws, whilst in reality our value 
of 7 is allowed to vary in the fits. However, we show in Table IbTI 
that using a fixed 7 = 1.8 only has a small affect on our ro values. 

One clear observation is that in both the volume-limited and 


magnitude-limited lo west ma s s sam ple the measured amplitude is 
much lower than the iLi et alj ( l2006h relation. As previously men¬ 
tioned, GAMA is unusually underdense at low redshifts, and the 
jack-knife error bars can considerably underestimate the true uncer¬ 
tainty at low redshift. The very low amplitude in the least massive 
samples is therefore likely sample variance. 

It appears at lower masses the dependence on ro with mass is 
stronger in our data than the fitting formulae. Some of this may be 
sample variance, but note this strong trend is not seen when only 
considering volume-limited samples. Indeed the volume-limited 
samples follow the expected dependence of clustering amplitude 
with mass. Another likely selection effect is that there is more 
redshift evolution in the magnitude-limited samples than expected 
from the passive evolution model. In Fig. [T^ we help to mitigate 
the effect of the GAMA r-band magnitude selection on our results 
by plotting ro against the median mass of the sample. Nonetheless, 
higher redshift samples will be biased toward different galaxies. 
As mentioned, more magnitude-limited sa mples will have lowe r 
clustering than volume-limited samples (e.g. iMeneux et al.ll2008h . 
As such the redshift evolution in clustering amplitude observed 
in the less massive galaxies may be down to these observational 
effects. Tests we conducted comparing clustering measurements 
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Figure 9. The projected two-point correlation function of real GAMA (black), G14 model (red) and L14 model (blue) galaxies as a function of stellar mass 
for different approximately volume-limited redshift slices. The real GAMA samples consist of at least 98% volume-limited galaxies, the mocks 89%. The 
measurements are divided by the reference power law defined in Section [T2l Error bars are from jack-knife resamplings for the real GAMA data and the 
scatter between the mock realisations for the mocks. Marked with P < 2% are samples for which the hypothesis that the data is a realisation of the model has 
less than a 2% probability, colour coded according to the model. 


from r < 19.8 and r < 21 versions of the mock catalogues support 
this idea (these samples have a much larger fraction of magnitude- 
limited galaxies than the samples where we found no significant 
differences between the original and deeper catalogues). 

At higher masses our measurements are not precise enough 
to determine if the galaxy clustering is evolving passively or not. 
It is likely that mass samples in these data can he used to probe 
redshift evolution, but this may be better done using more sophis¬ 
ticated modeling such as HOD fitting. We leave this to the HOD 
fitting of Palamara et al (in prep), while we focus on comparison 
to the models which have selection effects included via the use of 
lightcone. 

Returning to the model comparisons then, we see the ro val¬ 
ues of each model are very similar to one another, except in the high 
redshift slice where the L14 model has stronger clustering than the 
G14 model. The trends seen in the data are also reproduced by both 
models. In general the increase of ro with magnitude is slightly 
steeper in the model than the data. The steep increase of clustering 
with mass, at low masses, is expected to be related to the GAMA 
selection function. As this trend is qualitatively reproduced by the 
models, the models have some success in assigning the correct lu¬ 
minosities or colours to the sample galaxies. 


4.2.3 Summary of model comparisons for mass samples 


To summarize, the most obvious problems with the model predic¬ 
tions as a function of stellar mass are in the one-halo term regime, 
this points towards the physics of s atellite galaxies being a weak¬ 
ness in the model. As mentioned, in lCampbell et alJ 1 120151) instead 
of using the true model masses, they estimate the model masses 
from the predicted broad-band photometry. They find that this does 
affect the clustering as a function of mass. The effect is much 
stronger for the L14 model than the G14 model, and brings the 
small scale clustering into better agreement with the data. They 
still find that even by estimating the masses in this way, the models 
predict too much small-scale clustering. Also note that the typi¬ 
cally redder and fainter satellite galaxies are more sensitive to the 
r-band apparent magnitude cut (i.e. selection effects), so it may 
also be related to the model assigning wrong colours or luminosi¬ 
ties to satellite galaxies. This is also a possible explanation for the 
steeper increase in ro with mass compared to real data. 
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Figure 10. The ro of our power-law fits to the range 0.2 Mpc/i“^ < fp < 9.0 as a function of the median stellar mass of the sample divided 

by the characteristic mass of the stellar mass function (at 2 = 0). The errors come from fitting to the multiple jack-knife regions, or from the scatter from 
fitting to multiple mocks. Different panels give the GAMA measurements, the L14 model and the G14 model (left to right). Points have been colour coded b y 
median redshift, with samples defined by the same redshift cuts connected by lines. The dashed coloured lines gives the bias fitting formula of iLi et all i20Q6h. 
conver ted from bias to rg assuming a power law correlation function with 7 = 1.8 and evolved to different redshifts using the passive evolution model o dpryl 
<199d) . more details are given in the text. More massive galaxies are more clustered, a trend qualitatively reproduced by the models. Volume-limited samples 
are marked with a star, magnitude-limited samples are marked with a triangle. Note the magnitude-limited samples are likely to give lower values of ro than 
the corresponding volume-limited sample. As such the strong trend with redshift for the low mass samples should not be interpreted as clear evidence of real 
redshift evolution. This selection effect is accounted for when constructing the model catalogues, so the models and data can be fairly compared. 


4.3 Luminosity dependent clustering 

4.3.1 The full shape of the correlation functions 

In Fig.[TT]we show the clustering of galaxies as a function of lumi¬ 
nosity and redshift, in Fig. [T^ we show the corresponding volume- 
limited samples. For all of the redshift intervals we notice segre¬ 
gation between the faint and bright samples, with brighter galax¬ 
ies being more clustered. This trend is not seen in the mocks, in¬ 
deed the clustering of the L14 mock decreases slightly between 
—21.0 < Mr^h < —20.0 and —22.0 < Mr^h < —21.0. Because 
of this lack of dependence on luminosity, the mocks brighter than 
our sample —21.0 < Mr,h < —20.0 (which contains M* ff) all 
disagree very significantly with the data. For these bright galaxies 
the G14 model is closer to the data than the L14 model, but still 
very discrepant. 

For bins fainter than M* we see good agreement for many 
of the L14 model samples, in particular for samples fainter than 
Mr.h = —20. In contrast, some of the G14 model predictions 
for fainter luminosity samples show too much clustering on small 
scales (less than ~ 3 /i“^Mpc). At 0.24 < z < 0.35 and —21.0 < 
Mr.h < —20.0 the L14 model predicts too little small-scale clus¬ 
tering. Additionally, the most luminous mock samples under pre¬ 
dict the clustering, particularly on small scales (< 3/i“^ Mpc) at 
2 > 0.14. Apart from this most luminous sample, the small scale 
clustering predictions of the mocks are more successful here than 
in the mass samples. The volume-limited samples in Fig. [12] all 
show reasonable agreement, given the errors, with the model except 
the most luminous sample. In this sample, as with the magnitude- 
limited sample around this redshift, the models underpredict the 
clustering. 


4.3.2 Results from power law fits 

In Fig. [T^we plot the ro of our luminosity samples, as a function 
of the ratio of the m edian sample luminosity to = —20.6 

jLovedav et alJ[2oT^ . Taking this ratio is useful as it removes a 
dependence on how magnitudes are converted to luminosities. We 
see less variation in the clustering properties over this range than 
we did with mass. Volume-limited samples are plotted with a star, 
magnitude- limited samples are plotted with a triangle. In Fig. [T0| 
we plot the lZehavi et alj (1201 ih bias versus luminosity fitting for¬ 
mula as a dashed line, colour coded according to what redshift 
the relation has been passively evolved. Following Section [4.2.21 
to determine the rn .dm needed to convert galaxy biases from the 
IZehavi et alj ( 1201 ih formula to rp values, we use the quoted ro 
value from the IZehavi et alj ( 1201 ih power law fit to their sample 
containing L* as a reference. 

Before interpreting these data, we must once more consider 
the selection effects of the magnitude-limited samples. For these 
samples the effect is smaller, as the luminosity range is limited 
in the sample. Nonetheless, because of the colour-dependent k- 
corrections the colours of the samples can change due to the magni¬ 
tude selection. As before the selection effects may act to artificially 
enhance any redshift or luminosity trends, by making the clustering 
of fainter and more distant galaxies appear weaker than it would in 
a volume-limited sample. This certainly seems to be the case for 
the most luminous sample in the highest redshift range, where the 
volume-limited sample has a much larger amplitude than the mag¬ 
nitude limited one at the same redshift. Note the redshift range of 
the volume limited sample is very small, and as such it is likely 
to suffer from sample variance. The ro values of the other volume- 
limited samples seems to roughly agree with the magnitude-limited 
sample, but note this is for the subset of samples where producing 
a volume-limited version was possible. It was not possible to con¬ 
struct volume-limited samples for samples with very high fractions 
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Figure 11. The projected two-point correlation function of real (black) and G14 model (red) and L14 model (blue) galaxies as a function of redshift ( different 
columns) and luminosity (different rows). The measurements are divided by the reference power law defined in Section [52] Error bars are from jack-knife 
resamplings for the GAMA data, and from the scatter between realisations in the mock. Marked with P < 2% are samples for which the hypothesis that 
the data is a realisation of the model has less than a 2% probability, colour coded according to the model. The clustering of faint galaxies is successfully 
reproduced in the L14 model for the majority of samples. However, the models do not show as much clustering evolution with luminosity as the data, and are 
particularly discrepant for the brightest sample and on small scales. 


of magnitude-limited members, where selection effects are likely 
to be largest. 

The data show an increase in clustering with luminosity in all 
of the redshift slices. This is seen in many of the volume-limited 
samples, and so cannot solely be down to selection effects. Un¬ 
fortunately, as with the mass samples, this trend cannot be con¬ 
firmed in the two highest redshift slices with volume-limited sam¬ 
ples. Except for the lowest-luminosity sample our measure ments 
from real GAMA data agree with the IZehavi et alj ( 1201 ll) rela¬ 
tion in the lowest redshift slice, which is close to the redshift at 
which the relation was m easured. Our lowest luminosity measure¬ 
ment falls well below the lZehavi et alj j201 ih relation. As with the 
lowest- mass sample, becau se we know GAMA is under-dense at 
z < 0.1 ( [Driver et al J2011|) (which is above the median redshift of 
the faintest sample), a far more likely explanation is sample vari¬ 
ance unaccounted for by the jack-knife errors. 

Except for the brightest sample, the ro values seem to roughly 
follow the expectations of passive evolution. Of course, however, 
one has to worry about the selection effects. Using the volume- 
limited samples alone only gives some weak evidence of redshift 
evolution in ro between the lowest and second lowest redshift 
slices. 


The highest luminosity sample shows some evidence of not 
following the passive evolution model. The volume-limited sam¬ 
ple in particular shows stronger than expected clustering from a 
passive evolution model. However, given the size of the uncertain¬ 
ties the measurements are only I-2a away from the passive evolu¬ 
tion model. A lack of evolution, compared to that expected from 
passive evolution, would suggest fast bias evolution. This would 
be consistent with many other observations of more massive and 
brighter g alaxies showing li t tle evolution in clustering with red¬ 
shift (e.g. White et alj |2Q0^ : [Brown et alj [200 A ICoil et ahlbOOSl : 
iMeneux et al.ll2008h. 


Moving back to our focus of model comparisons, the mod¬ 
els show little dependence on clustering amplitude and luminos¬ 
ity. Both models predict a similar trend of ro versus luminosity. 
This lack of l umino sity dependence is also in agreement with the 
IZehavi et alj (1201 ll) bias fitting formula, up to the brightest sam¬ 
ple where the L14 model ro values falls below the fitting formula. 
The models both predict much smaller ro values than the data in 
the highest luminosity bin, the models, unlike the data, also show 
evidence of evolution in the clustering with redshift in this sample. 
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Figure 12. The projected two-point correlation function of real (black) and G14 model (red) and L14 model (blue) galaxies for different volume-limited 
luminosity samples. The measurements are divided by the reference power law defined in Section [21 Error bars are from jack-knife resamplings for the 
GAMA data, and from the scatter between realisations in the mock. Marked with P < 2% are samples for which the hypothesis that the data is a realisation 
of the model has less than a 2% probability, colour coded according to the model. 


4.3.3 Summary of model comparisons for luminosity samples 


The models do a reasonable job at reproducing the amplitude of 
clustering as a function of luminosity for the samples fainter than 
L*, and the L14 model is a good fit to many of the fainter sam¬ 
ples. However, the models fail to reproduce the amplitude of our 
brightest sample. This under prediction particu larly affects small- 
scale clustering. Also, recall Ixim et aP ( l2009l) found the depen- 
dence on luminosity was not correct for the Bower et"^ l l2006l) 
model. Note that semi-analytic models have been shown to have a 
dependence on clustering with luminosity, but this evolution only 
star ts to become apparen t at the highest luminosity probed here (see 
e.g. lNorberg et al.ll200lh . Again, the disagreements are particularly 
large at small scales, suggesting the modelling of satellite galaxies 
needs to be improved in the models. 


4.4 Colour dependent clustering 

Fig-[M]shows the clustering of our samples of red and blue GAMA 
galaxies (in black). In all redshift intervals we probe, red galaxies 
are more clustered than blue galaxies. We also note that the red 
galaxy correlation functions have their strongest clustering, rela¬ 
tive to the reference line, at scales less than around 2 h~^ Mpc. 
This can be interpreted as red galaxies predominantly being in 


larger haloes with more satellite galaxies. The blue galaxy sam¬ 
ples also show a relative increase in clustering on very small scales 
(< 0.3 h~^ Mpc), but note that their clustering is still weaker than 
the reference power law. 

In each panel of Fig. M we indicate the absolute magnitude 
range enclosing the central 68% of the data. From Fig. m we can 
draw the following natural conclusions. Firstly, for a given redshift 
slice, red galaxies are typically brighter than blue ones. Addition¬ 
ally, the absolute magnitude range of the central 68-percentile re¬ 
duces with redshift and gets brighter with redshift. Finally, at fixed 
redshift, the change in clustering between red and blue galaxies 
is significantly larger than can be explained by any luminosity (or 
stellar mass) dependence of clustering alone. The models also show 
red galaxy samples have a larger amplitude and steeper 2PCFs than 
blue galaxies. At low redshifts (z < 0.24) the L14 model agrees 
with the measurements whilst the G14 model has too much small- 
scale clustering. In contrast the G14 model shows better agreement 
in the two high redshift slices, whilst the L14 model now has too 
little small-scale clustering. 

The blue mock galaxy samples all have too low an amplitude, 
particularly at small scales. The size of the values can leave us in 
no doubt that neither model successfully reproduces the clustering 
of blue galaxies. As the small scale measurements are particularly 
discrepant, it is possible that too few satellite galaxies are blue. 
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Figure 13. The tq of our power law fits to the range 0.2 < rp < 9.0/i~^Mpc as a function of the median luminosity of the sample, divided by the 
characteristic luminosity of the data at 2 : = 0. Pure luminosity evolution has been accounted for with Q=1.45. The error bars give the error from fitting to the 
multiple jack-knife regions. Different panels give the GAMA measurements, the G14 model and the L14 model (left to right). Points have been colour coded 
by median redshift, with samples defined by the same redshift cuts con nected by lines. Volu me-limited samples are plotted as stars, magnitude-limited samples 
as triangles. The dashed coloured lines gives the bias fitting formula of lZehavi~ with fjQ.s = 0. 8, converte d from bias to 7*0 assuming a power law 

correlation function with 7 = 1.8 and evolved to different redshifts using the passive ev olution model 0flR^h99ell . more details are given in the text. The 
IZehavi et al.N20l'l|) fitting formula is only derived from the range 0.16 < Lr/L* < 6.3 iZehavi et al.l201 ilk but the disagreement at low luminosity is more 
likely to be down to sample variance underestimated in the jack-knife en'ors. Brighter galaxies are more clustered in the data, a trend not reproduced in the 
models. Note the magnitude-limited samples are likely to be underestimates of the true vq. As such the trends with redshift for the low luminosity samples 
should be interpreted with caution (see the text). In the highest luminosity samples, three volume-limited samples demonstrate that ro is evolving slower than 
the passive evolution model (i.e. the bias is evolving faster). 


As these samples are magnitude-limited, another possibility is the 
model predicts blue satellites that are too faint. The predictions for 
blue galaxies being incorrect may not be too surprising, when one 
considers that the models predict a clear sequence of blue galaxies 
not seen in the data (Fig. n. The models also do not reproduce 
the increase in clustering relative to the power law seen at very 
small scales (< 0.3 h~^ Mpc). This could relate to an increase 
in star formation rate for close pairs of galaxies. This possibility 
will be explored in a future paper using star-formation-rate-selected 
samples of GAMA galaxies (Gunawardhana et al. in prep). 


5 DISCUSSION AND CONCLUSIONS 

We have studied the projected two-point correlation function of 
galaxies in GAMA. To do this we used a modified version of the 
IColel ( I2OI ll) approach to generate random catalogues; this method 
resulted in a set o f rand o m poi nts with all of the properties of the 
real galaxies. The IColel ( 1201 ih approach allows sample selection 
cuts to be applied to both the data and the random catalogue, allow¬ 
ing the measurement of galaxy clustering as a function of diverse 
galaxy properties. Our modification, the inclusion of a window to 
limit the difference between the initial and cloned redshift, has the 
potential to limit the effects of galaxy evolution on random cata¬ 
logues. This method should be followed up in later work, to test 
how effective it is on a wider selection of galaxy samples (e.g. not 
just selected on optical photometry) and to optimize the size of win¬ 
dow used. 

We compared volume-limited samples from SDSS of 
IZehavi et al.l ( 12011[) to volume-limited samples in GAMA, both 
with the redshift cuts of SDSS and with redshift cuts appropri¬ 
ate for the deeper, GAMA magnitude limit. We find good agree¬ 


ment with SDSS for our samples using the SDSS redshift limits. In 
our —23.0 < Mjl < —22.0 sample and in our faintest sample, 

— 19.0 < < —18.0, both with the GAMA redshift limits, 

we find some disagreement between SDSS and GAMA. Given the 
expected uncertainties in our jack-knife covariance matrices, this 
could indicate the errors were underestimated. Our faintest sample, 

— 19.0 < < —18.0, shows less cluste ring than the SDSS 

data. Particularly relevant for this sample is the [Driver et ahl ( 1201 ih 
observation that GAMA is under-dense at z < 0.1, which would 
indeed lead to lower clustering. 


We have observed that more luminous, more massive and red¬ 
der galaxies are more strongly clustered, in redshift slices between 
z = 0 and z = 0.5. Though in the two highest redshift slices 
the GAMA selection function complicates the interpretation of this 
trend. More massive and luminous galaxies being more clustered 
is in agreement with previous measurements th at show these trends 
exist at lower and higher r edshift ranges (e.g. IZehavi et al. 


Christodoulou et aL| 201^ UetaU [201 2l : Ide la Torre et al. 


2011 


2013 


Marulli et alj2013l ; lGuo et al.l2014l) . and in a broad redshift bin en¬ 


compassing all of our slices dSkibba et alj|2014h . We also find that 
red galaxies have steeper correlation functions than blue galaxies 
in these redshift slic es, again in agreement with clus tering mea- 
sure d at higher (e.g. l(3oil et alj|2008l ; IGuo et alj|2014l) and lower 
(e.g. lZehavi et al.ll2011 ) redshifts. 


We fit power laws to our measurements and see an evolution in 
the apparent clustering strength of galaxies with redshift for sam¬ 
ples less massive than MQh~^. We also evidence of this oc¬ 

curring for galaxies less luminous than around L*. The evolution 
is in the direction of higher redshift galaxies being less clustered. 
Note however, many of our samples are not volume limited and so 
will be affected by selection effects that are very likely to mimic 
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Figure 14. The projected correlation function of red galaxies (top) and blue galaxies (bottom) in different redshift slices. Also shown is the clustering of mock 
galaxies for the L14 model (blue) and the G14 model (red). The results have been divided by the reference power law defined in SectionMarked with 
P < 2% are samples for which the hypothesis that the data is a realisation of the model has less than a 2% probability, colour coded according to the model. 


this evolution. The volume-limited samples alone do show some 
evidence of higher redshift samples of the same luminosity or mass 
being less clustered, but this evidence is rather weak considering 
the size of the uncertainties. 

In our most luminous sample (~ 3L*), which is volume- 
limited in three of the redshift slices, we see little evolution in clus¬ 
tering amplitude between samples with median redshifts of 2 ~ 0.1 
and 2 ~ 0.4. Given our uncertainties, there is some weak evidence 
that this is less evolution than that expected from a simple, passive 
evolution model. 

A lack of evolution in the clustering amplitude of brighter and 
more massive galaxies has been observed by other authors, and 
taken as evidence of a fast evolution of the bias of these object s 
dWhite et alj|2007t iBrown et alJboOSi: iMeneux et al1l2008l. l2009t) . 

This fast bias evolut ion has been connected to the merging and dis¬ 
ruption of satellites dWhite et al.ll2007l : [Brown et alJl2008 ). 

Whilst the interpretation of these measurements is compli¬ 
cated by selection effects in the magnitude-limited samples, the 
strength and focus of this work is model comparisons. We use 
a model that utilizes a lightcone, and assign fc-corrections to the 
model galaxies in such a way to mimic the selections function of 
the survey. We find that the L14 and G14 semi-analytic models 
successfully predict the trends of clustering as a function of stel¬ 
lar mass. Both models have similar predictions for this trend, sug¬ 
gesting it is insensitive to differences in the adopted physics. In 


detail there are places where the models disagree with the data, 
often in the regime of the one-halo term, which is too strong. 
ICampbell et^ (12015h notes that estimating model masses from the 
colours, rather than taking them from the model, helps alleviate this 
problem for the L14 model. 

We find the increase of clustering seen between our samples 
—21.0 < Mr^h < —20.0 and —22.0 < Mr^h < —21.0 is not re¬ 
produced by the models. This discrepancy is present in both mod¬ 
els, but is a bigger problem for the L14 model. For the highest lumi¬ 
nosity sample (~ 3L*) the models also predict too little small-scale 
clustering. The clustering of fainter galaxies is better modelled in 
the data, with the L14 model being particularly successful at repro¬ 
ducing the clustering of the Mr,h > —20 galaxy samples. 

We have also shown that the trends of redshift evolution with 
mass and luminosity are qualitatively reproduced by the models 
(Fig. [To] and Fig. [Tsl l. in that the more massive and more lumi¬ 
nous samples show less redshift evolution. The redshift trends are 
most likely successfully reproduced thanks to the modelling of 
the GAMA apparent magnitude limit, showing the importance of 
applying observational effects when comparing galaxy formation 
models to data. 

The clustering of red galaxies on large scales is successfully 
reproduced by the models. At 2 > 0.35 the L14 model under pre¬ 
dicts the clustering on small scales and at 2 < 0.35 the G14 model 
over predicts clustering on small scales. Both models under-predict 
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the clustering of blue galaxies, particularly on scales less than a 
few Mpc. The colour-magnitude diagram of the models is also very 
different to that of the data, the models have a much more clearly 
bimodal colour distribution. 

These models give different predictions from one another for 
samples defined by mass, photometric properties and colour. This 
highlights the importance of testing the clustering of models using 
a varie ty of properties for sample selection, or, as ICampbell et alj 
( l2015h suggests, inferring the mass of the model galaxies from their 
photometric properties. 

The models struggle to reproduce the clustering at small 
scales. Therefore, a suggested route to improvement is modifying 
the physics affecting satellite galaxies. Several authors have ex¬ 
perimented with adding processes which remo ve satellite galax¬ 
ies, or modify their luminosity or colour (e.g Font et all 20081 : 

' l l2009l) addition¬ 


les, or modity their luminosity or c olon 
iKim et ani2009l : IContreras et alj[20l^ . Ie 


im et al 


ally found that processes affecting satellite galaxies can also change 
the dependence of clustering on luminosity. Modifying the satel¬ 
lite galaxy physics could also have an impact of the colour dis- 
tributio n and improve the clustering predictions of blue galaxies. 
Indeed, Ipont et aP bOOSli found bluer satellites when applying a 
more realistic model of the stripping of gas from infalling satellite 
galaxies. Modifications to the satellite galaxy physics could also af- 
fect the over abundance of faint red galaxies in the models noted by 
iMcNaught-Roberts et al.l ( l2014h . or the excess of high-m ultiplicity 
groups of galaxies observed by iRobotham et al.l ( l201lh . Finally, 
changing the colours of satellite galaxies can have knock-on ef¬ 
fects for magnitude-limited samples, as model galaxies with red 
SEDs are fc-corrected out of these samples at lower redshifts than 
blue galaxies. Such modifications should also be tested to see if 
they reproduce the lack of redshift evolution in the brightest sam¬ 
ples where the physics affecting satellites galax ies is expected to 
play a role dWhite et al.ll2007l : [Brown et alfcOOSh . 

In conclusion, the G14 and L14 models reproduce the trends 
of clustering with mass and successfully predict the amplitude of 
clustering for galaxies fainter than L*. They do, however, need im¬ 
provement at small scales - in particular blue galaxies need stronger 
small-scale clustering. They also need to be improved in order to re¬ 
produce the clustering of galaxies slightly brighter than L*, i.e. for 
our samples with median r-band magnitudes around Mr^h ~ 21.3. 
When these measurements are compared to models one must be 
careful to account for the selection effects in the magnitude-limited 
samples, for a realistic mock catalogue this could be as simple as 
applying an r < 19.8 apparent magnitude cut. An alternative is 
to choose the samples which are volume-limited. Indeed, in a fu¬ 
ture GAMA paper we intend to use the methodology here to pro¬ 
duce volume-limited magnitude and mass threshold selected sam¬ 
ples and carry out a HOD analysis on them (Palamara in prep.). In 
order to facilitate further model tests against these measurements, 
our results will be made available online at publicatior0. 
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APPENDIX A: THE EFFECT OF P, Q AND WINDOWED 
RANDOMS ON THE CORRELATION FUNCTION 

In this paper we utilize a random catalogue generated using a new, 
windowed approach. We also use a different combination of the 
luminosity function evolution parameters Q — 1.45, P — 0.0 
to a different set of parameters taken from an earlier draft of 
iLovedav et al.l (l2015h : Q — 0.81, P = 1.45. To test for any biases 
introduced by this, we compare results of using three different ran¬ 
dom catalogues: the randoms used in our measurements, the same 
P and Q as our random catalogue but without the window function 
and a different P and Q with and without the window function. For 
this purpose we use a mass selected sample, as this avoids changes 
in the sample definition arising due to a different value of Q. We 
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Figure Al. The projected 2PCF of 10.00 < log;LQ < 10.50 

galaxies in two different redshift slices, as indicated. The different lines 
represent different methods of generating a random catalogue with which 
to measure clustering, each is labelled in the legend. It can be seen that the 
our choice of random catalogue does not have a significant effect on our 
measurements. 
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choose the 10.00 < logj^Q < 10.50 mass bin, and use 

the redshift ranges 0.14 < 2; < 0.24 and 0.24 < 2 < 0.35. These 
redshift slices encompass a large under-density in the n{z), where 
one might expect to see the largest differences between the random 
catalogues. 

In Fig. lAll we show the results. The most difference appears at 
large-scales; with the measurements using the windowed catalogue 
showing a lower clustering amplitude than the results from the two 
random catalogues without windowing. The two unwindowed re¬ 
sults also show some deviation on large-scales. Note, however, that 
the differences are much smaller than the error bars for both of the 
samples, so much so that the conclusions of this paper will not be 
affected. 


APPENDIX B: POWER-LAW LIT TABLES AND MOCK 
SAMPLE PROPERTIES 

In Table IBTI we show the results from our power law fits to our 
samples of GAMA galaxies. Uncertainties are given in brackets, 
and were derived from fitting power laws to each jack-knife real¬ 
isation separately. Also given are the ro values found from fixing 
the slope of the power law, 7, rather than allow it to vary freely. 
The difference in the recovered ro values is small. 

Tables lB2l and lB3| give the sample properties, and standard de¬ 
viations (in brackets) of one realisation of the G14 and L14 models 
respectively. The magnitudes have been adjusted in order to match 
their luminosity functions to the real GAMA data (see l2.3.3l l. as 
such the sample sizes are much more similar to the real data for 
luminosity-defined samples than mass- or colour-defined samples. 
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Table Bl. The recovered values from power law fits to the clustering, ^0,7 fixed fhe results of fitting the data with a fixed value of 7 = 1 . 8 . Jack-knife 
errors are given in brackets. We see red galaxies have steeper 2PCFs, and that the amplitude of the 2PCF increases with magnitude and mass, regardless as 
to whether 7 is also fit or not. Samples marked with asterisks are magnitude limited, so should be treated with careful consideration of the GAMA selection 
function. Samples without asterisks are volume limited, and can be treated at face value. 


Sample 

'^min 

^max 

ro 

7 

^0,7 fixed 

-18.00 < < -17.00 

0.02 

0.14* 

3.68 (0.31) 

1.81 (0.09) 

3.69 (0.27) 

-19.00 < Mr^H < -18.00 


0.14* 

4.57 (0.37) 

1.74 (0.06) 

4.44 (0.26) 



0.24* 

4.51 (0.69) 

1.84 (0.05) 

4.52 (0.70) 

-20.00 < < -19.00 

0.02 

0.14 

5.35 (0.51) 

1.70 (0.07) 

5.07 (0.36) 


0.14 

0.24* 

4.91 (0.22) 

1.78 (0.03) 

4.86 (0.17) 


0.24 

0.35* 

4.05 (0.17) 

1.69 (0.05) 

4.09 (0.18) 

-21.00 < Mr,h < -20.00 

0.02 

0.14 

5.58 (0.75) 

1.76 (0.07) 

5.45 (0.53) 


0.14 

0.24 

5.39 (0.28) 

1.78 (0.04) 

5.33 (0.19) 


0.24 

0.35* 

5.04 (0.09) 

1.76 (0.02) 

4.98 (0.07) 


0.35 

0.50* 

4.82 (0.27) 

1.62 (0.06) 

4.84 (0.26) 

-22.00 < Mr^h < -21.00 


0.14 

6.59 (0.58) 

1.77 (0.08) 

6.49 (0.44) 


0.14 

0.24 

6.12(0.32) 

1.87 (0.05) 

6.30 (0.25) 


0.24 

0.35 

6.47 (0.18) 

1.79 (0.02) 

6.45 (0.16) 


0.35 

0.50* 

6.38 (0.25) 

1.81 (0.03) 

6.40 (0.20) 

-23.00 < MO I < -22.00 

0.01 

0.50 

9.39 (0.50) 

1.99 (0.10) 

9.56 (0.44) 

-22.00 < MO I < -21.00 

0.01 

0.38 

6.65 (0.10) 

1.77 (0.01) 

6.60 (0.09) 

-21.00 < < -20.00 

0.01 

0.26 

5.72 (0.28) 

1.75 (0.03) 

5.55 (0.20) 

-20.00 < AfO i < -19.00 

0.01 

0.18 

5.27 (0.42) 

1.74 (0.04) 

5.08 (0.31) 

-19.00 < M7 < -18.00 

0.01 

0.12 

3.97 (0.33) 

1.87 (0.06) 

4.09 (0.27) 

8.50 < logj^Q Mf /MQh~^ < 9.50 

0.02 

0.05 

3.28 (0.42) 

2.24 (0.21) 

4.86 (0.38) 


0.02 

0.14* 

4.00 (0.22) 

1.73 (0.05) 

3.95 (0.21) 


0.14 

0.24* 

3.70 (0.19) 

1.66(0.03) 

3.65 (0.17) 

9.50 < logioM*/M©/i-2 < 10.00 

0.02 

0.14 

5.84 (0.64) 

1.77 (0.07) 

5.72 (0.37) 


0.14 

0.24* 

4.81 (0.24) 

1.77 (0.03) 

4.75 (0.20) 


0.24 

0.35* 

3.44 (0.11) 

1.64(0.05) 

3.64 (0.09) 

10.00 < logio M*/M0h-2 < 10.50 

0.02 

0.14 

5.93 (0.71) 

1.82(0.07) 

6.01 (0.54) 


0.14 

0.18 

5.94 (0.46) 

1.86 (0.05) 

6.35 (0.36) 


0.14 

0.24* 

5.67 (0.27) 

1.85 (0.04) 

5.86 (0.17) 


0.24 

0.35* 

4.81 (0.13) 

1.73 (0.02) 

4.71 (0.11) 


0.35 

0.50* 

4.49 (0.20) 

1.69(0.04) 

4.65 (0.19) 

10.50 < logio M*/M©h-2 < 11.00 

0.02 

0.14 

6.84 (0.62) 

1.75 (0.09) 

6.69 (0.47) 


0.14 

0.24 

6.37 (0.34) 

1.93 (0.04) 

6.78 (0.30) 


0.24 

0.29 

5.76 (0.28) 

1.93 (0.03) 

6.65 (0.24) 


0.24 

0.35* 

6.37 (0.20) 

1.84 (0.02) 

6.50 (0.17) 


0.35 

0.50* 

6.15 (0.22) 

1.78(0.03) 

6.10(0.16) 

11.00 < logio M*/M©h-2 < 11.50 


0.35 

8.42 (0.37) 

1.76(0.03) 

8.33 (0.34) 


0.35 

0.37 

8.06(1.33) 

2.14(0.11) 

10.88(1.19) 



0.50* 

8.17 (0.47) 

1.88(0.05) 

8.38 (0.44) 

Red 


0.14* 

6.25 (0.53) 

2.03 (0.07) 

7.13 (0.45) 


0.14 

0.24* 

6.37 (0.31) 

1.96 (0.03) 

7.04 (0.25) 


0.24 

0.35* 

6.49(0.15) 

1.90(0.01) 

6.79 (0.14) 


0.35 

0.50* 

7.50 (0.29) 

1.89 (0.02) 

7.82 (0.25) 


0.02 

0.14* 

3.34 (0.14) 

1.52 (0.02) 

3.03 (0.10) 

0.14 

0.24* 

3.89 (0.14) 

1.63 (0.03) 

3.63 (0.10) 

0.24 

0.35* 

4.09(0.11) 

1.64(0.02) 

4.00 (0.10) 

0.35 

0.50* 

4.80 (0.10) 

1.69 (0.03) 

4.85 (0.10) 


Blue 
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Table B2. Different mock galaxy samples, sample sizes and median properties for one G14 model realisation of GAMA. The subscript ‘med’ indicates the 
values are medians. Where maximum redshift has an asterisk, the sample is not volume limited. In samples without an asterisk at least 95% of the members 
are volume limited. Values in brackets are rms scatter. 



Sample 

'^min 

■^max 

-^gals 

•^med 


loglo(^«/'^0^ ^)med 

(9 ^)o,med 


-18.00 < Mr,h < -17.00 



2965 


-17.46 (0.28) 

8.83 (0.26) 

0.29 (0.17) 


on 

0.14* 


la 

-17.59 (0.27) 

8.83 (0.25) 

0.27(0.17) 


-19.00 < Mr,h < -18.00 

0.02 

0.11 

7402 

0.08 

-18.48 (0.29) 

9.20 (0.30) 

0.30 (0.18) 


0.02 

0.14* 

12979 

0.10 

-18.55 (0.27) 

9.21 (0.30) 

0.29(0.18) 



0.14 

0.24* 

2894 

0.15 

-18.85 (0.12) 

9.28 (0.33) 

0.30(0.19) 


-20.00 < Mr^h < -19.00 

0.02 


11689 

0.11 

-19.47 (0.29) 

9.71 (0.37) 

0.43 (0.20) 



0.14 


7029 

0.15 

-19.45 (0.28) 

9.75 (0.39) 

0.51 (0.20) 



0.14 


26387 

0.19 

-19.61 (0.27) 

9.78 (0.39) 

0.49 (0.21) 



0.24 


2343 

0.25 

-19.91 (0.09) 

9.65 (0.29) 

0.27(0.17) 


-21.00 < Mr^h < -20.00 


0.14 

6968 

0.11 

-20.40 (0.28) 

10.26 (0.36) 

0.61 (0.21) 


0.14 

0.24 

22285 

0.20 

-20.39 (0.28) 

10.28 (0.35) 

0.61 (0.21) 



0.24 

0.35* 

35746 

0.28 

-20.56 (0.26) 

10.24 (0.36) 

0.50(0.21) 



0.35 

0.50* 

3644 

0.37 

-20.87 (0.14) 

10.05 (0.28) 

0.27(0.17) 


-22.00 < < -21.00 


0.14 

1835 

0.11 

-21.28 (0.25) 

10.72 (0.32) 

0.70 (0.20) 


0.14 

0.24 

5838 

0.20 

-21.28 (0.25) 

10.73 (0.32) 

0.70 (0.20) 



0.24 

0.35 

15354 

0.30 

-21.29 (0.25) 

10.72 (0.32) 

0.66 (0.20) 



0.35 

0.37 

3316 

0.36 

-21.29 (0.26) 

10.72 (0.32) 

0.64 (0.20) 



0.35 

0.50* 

17073 

0.40 

-21.41 (0.26) 

10.64 (0.35) 

0.48 (0.21) 


8.50 < M* < 9.50 

0.02 

0.05 

3066 


-17.49 (0.89) 

8.92 (0.28) 

0.43 (0.16) 



0.02 

0.14* 

20046 


-18.34 (0.76) 

9.11 (0.26) 

0.27(0.14) 



0.14 

0.24* 

7525 

m 

-19.23 (0.30) 

9.37 (0.12) 

0.25 (0.07) 


9.50 < M, < 10.00 

0.02 


10713 

0.11 

-19.44(0.63) 

9.73 (0.14) 

0.54 (0.16) 



0.14 


19616 

0.19 

-19.79 (0.47) 

9.77 (0.14) 

0.35 (0.15) 



0.24 


12697 

0.28 

-20.34 (0.28) 

9.84 (0.12) 

0.27 (0.08) 


10.00 < M, < 10.50 

0.02 


5746 

0.11 

-20.18 (0.58) 

10.20 (0.14) 

0.66 (0.15) 



0.14 

0.18 

5155 

0.16 

-20.14(0.57) 

10.20 (0.14) 

0.66 (0.15) 



0.14 

■is [SI 

18085 

0.20 

-20.17(0.53) 

10.22 (0.14) 

0.65 (0.15) 



0.24 


19509 

0.28 

-20.72 (0.42) 

10.25 (0.15) 

0.44(0.16) 



0.35 

■ 1 

8752 

0.40 

-21.23 (0.30) 

10.24 (0.14) 

0.29 (0.09) 


10.50 < M, < 11.00 

0.02 

0.14 

3455 

0.11 

-20.71 (0.69) 

10.73 (0.14) 

0.76 (0.11) 



0.14 

0.24* 

11227 

0.20 

-20.69 (0.61) 

10.71 (0.14) 

0.75 (0.11) 



0.24 

0.29* 

8747 

0.27 

-20.83 (0.43) 

10.71 (0.14) 

0.72(0.11) 



0.24 

0.35* 

18649 

0.29 

-20.94 (0.40) 

10.72 (0.14) 

0.71 (0.11) 



0.35 

0.50* 

7885 

0.39 

-21.39 (0.33) 

10.80 (0.15) 

0.67(0.16) 


11.00 < M* < 11.50 


0.35* 

3114 


-21.75 (0.50) 

11.09(0.09) 

0.75 (0.05) 




0.37 

609 

^^B 

-21.79 (0.32) 

11.10(0.09) 

0.74 (0.05) 




0.50* 

4377 


-21.84 (0.28) 

11.11 (0.09) 

0.73 (0.05) 


Red 

0.02 

0.14* 

19247 

0.10 

-19.23 (1.37) 

9.87 (0.65) 

0.63 (0.10) 

(9- 

r)o + 0.03(Mr,h - M*,^) > 0.498 

0.14 

0.24* 

31606 

0.19 

-20.07 (0.70) 

10.35 (0.38) 

0.69 (0.09) 


0.24 

0.35* 

28164 

0.28 

-20.80 (0.51) 

10.67 (0.28) 

0.71 (0.08) 



0.35 

0.50* 

10034 

0.39 

-21.53 (0.38) 

10.97 (0.21) 

0.72 (0.06) 


Blue 

0.02 

0.14* 

24008 

0.09 

-18.72(1.27) 

9.18 (0.53) 

0.26 (0.07) 

(9 - 

r)o + 0.03(M^,^ - < 0.498 

0.14 

0.24* 

26135 

0.19 

-19.88(0.70) 

9.68 (0.34) 

0.30 (0.08) 


0.24 

0.35* 

26245 

0.29 

-20.63 (0.49) 

10.01 (0.28) 

0.30 (0.09) 



0.35 

0.50* 

12622 

0.40 

-21.25 (0.40) 

10.26 (0.26) 

0.29(0.10) 
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Table B3. Different mock galaxy samples, sample sizes and median properties for one L14 model realisation of GAMA. The subscript ‘med’ indicates the 
values are medians. Where maximum redshift has an asterisk, the sample is not volume limited. In samples without an asterisk at least 95% of the members 
are volume limited. Values in brackets are rms scatter. 



Sample 

'^min 

■^max 

-^gals 

•^med 


loglo(^«/'^0^ ^)med 

i.9 ^)o,med 


-18.00 < Mr,h < -17.00 



2982 


-17.46 (0.28) 

8.82 (0.30) 

0.41 (0.15) 


on 

0.14* 

5922 

la 

-17.59 (0.27) 

8.82 (0.29) 

0.39 (0.15) 


-19.00 < Mr,h < -18.00 

0.02 

0.11 

7348 

0.08 

-18.49 (0.29) 

9.21 (0.31) 

0.42 (0.15) 


0.02 

0.14* 

12815 

0.10 

-18.55 (0.27) 

9.22 (0.31) 

0.41 (0.15) 



0.14 

0.24* 

2887 

0.15 

-18.85 (0.12) 

9.27 (0.30) 

0.41 (0.16) 


-20.00 < Mr^h < -19.00 

0.02 


11665 

0.11 

-19.47 (0.29) 

9.65 (0.35) 

0.51 (0.16) 



0.14 


7242 

0.15 

-19.45 (0.29) 

9.67 (0.37) 

0.53 (0.17) 



0.14 


26443 

0.19 

-19.61 (0.27) 

9.70 (0.38) 

0.53 (0.17) 



0.24 


2560 

0.25 

-19.91 (0.09) 

9.62 (0.26) 

0.36(0.15) 


-21.00 < Mr^h < -20.00 


0.14 

6934 

0.11 

-20.41 (0.28) 

10.18 (0.40) 

0.62 (0.18) 


0.14 

0.24 

22310 

0.20 

-20.40 (0.28) 

10.22 (0.40) 

0.63 (0.18) 



0.24 

0.35* 

35739 

0.28 

-20.56 (0.27) 

10.16(0.41) 

0.56(0.19) 



0.35 

0.50* 

3563 

0.37 

-20.86 (0.15) 

10.00(0.31) 

0.37(0.16) 


-22.00 < < -21.00 


0.14 

1833 

0.11 

-21.28 (0.25) 

10.63 (0.39) 

0.65 (0.18) 


0.14 

0.24 

5843 

0.20 

-21.28 (0.25) 

10.65 (0.38) 

0.64(0.18) 



0.24 

0.35 

15349 

0.30 

-21.28 (0.25) 

10.63 (0.38) 

0.61 (0.18) 



0.35 

0.37 

3399 

0.36 

-21.29 (0.26) 

10.64 (0.38) 

0.59 (0.17) 



0.35 

0.50* 

17320 

0.40 

-21.42 (0.26) 

10.47 (0.39) 

0.49(0.18) 


8.50 < M* < 9.50 

0.02 

0.05* 

3077 


-17.46 (0.94) 

8.92 (0.28) 

0.55 (0.14) 



0.02 

0.14* 

19994 


-18.36 (0.80) 

9.12(0.26) 

0.38(0.12) 



0.14 

0.24* 

9080 

m 

-19.29 (0.31) 

9.36 (0.12) 

0.36 (0.08) 


9.50 < M, < 10.00 

0.02 

0.14* 

11290 

0.11 

-19.45 (0.71) 

9.72 (0.14) 

0.58 (0.13) 



0.14 

0.24* 

20525 

0.19 

-19.86 (0.50) 

9.76 (0.14) 

0.44(0.12) 



0.24 

0.35* 

15610 

0.28 

-20.39 (0.31) 

9.82 (0.12) 

0.38 (0.09) 


10.00 < M, < 10.50 

0.02 


5550 

0.11 

-20.09 (0.66) 

10.19(0.14) 

0.68 (0.11) 



0.14 

0.18 

5102 

0.16 

-20.06 (0.65) 

10.20 (0.14) 

0.68 (0.12) 



0.14 

■is [SI 

16659 

0.20 

-20.14(0.61) 

10.21 (0.14) 

0.67(0.12) 



0.24 


18009 

0.29 

-20.82 (0.47) 

10.22 (0.14) 

0.53 (0.12) 



0.35 

■ 1 

9983 

0.40 

-21.33 (0.33) 

10.22 (0.14) 

0.42 (0.09) 


10.50 < M, < 11.00 

0.02 

0.14 

2806 

0.11 

-20.71 (0.58) 

10.71 (0.14) 

0.78 (0.10) 



0.14 

0.24* 

9615 

0.20 

-20.66 (0.53) 

10.72 (0.14) 

0.78(0.10) 



0.24 

0.29* 

7770 

0.27 

-20.76 (0.44) 

10.72 (0.14) 

0.76(0.10) 



0.24 

0.35* 

16067 

0.29 

-20.88 (0.42) 

10.75 (0.14) 

0.75 (0.10) 



0.35 

0.50* 

6316 

0.39 

-21.47(0.38) 

10.78 (0.15) 

0.62(0.14) 


11.00 < M* < 11.50 


0.35* 

4280 


-21.42 (0.51) 

11.09(0.09) 

0.79 (0.07) 




0.37 

760 


-21.50 (0.35) 

11.10(0.09) 

0.78 (0.06) 




0.50* 

3912 


-21.68 (0.35) 

11.12(0.10) 

0.77 (0.08) 


Red 

0.02 

0.14* 

20986 

0.10 

-19.07(1.38) 

9.80 (0.65) 

0.65 (0.09) 

( 9 - 

r)o + 0.03(Mr,h - > 0.548 

0.14 

0.24* 

31557 

0.19 

-20.06 (0.67) 

10.29 (0.41) 

0.70 (0.09) 


0.24 

0.35* 

27728 

0.28 

-20.77 (0.47) 

10.68 (0.32) 

0.73 (0.08) 



0.35 

0.50* 

7944 

0.39 

-21.46 (0.36) 

10.98 (0.24) 

0.75 (0.07) 


Blue 

0.02 

0.14* 

22218 

0.10 

-18.80(1.30) 

9.18 (0.53) 

0.36 (0.07) 

(9 - 

r)o + 0.03(M^,^ - A4*^) < 0.548 

0.14 

0.24* 

26243 

0.19 

-19.88(0.73) 

9.62 (0.33) 

0.39 (0.08) 


0.24 

0.35* 

26870 

0.29 

-20.65 (0.54) 

9.95 (0.29) 

0.40 (0.09) 



0.35 

0.50* 

14824 

0.40 

-21.34 (0.43) 

10.22 (0.29) 

0.41 (0.10) 






















